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INTRODUCTION 


The Institute for Computer Applications in Science and Engineering (ICASE) is oper- 
ated at the Langley Research Center (LaRC) of NASA by the Universities Space Research 
Association (USRA) under a contract with the Center. USRA is a nonprofit consortium of 
major U. S. colleges and universities. 

The Institute conducts unclassified basic research in applied mathematics, numerical 
analysis, and computer science in order to extend and improve problem-solving capabilities 
in science and engineering, particularly in aeronautics and space. 

ICASE has a small permanent staff. Research is conducted primarily by visiting scientists 
from universities and from industry, who have resident appointments for limited periods of 
time, and by consultants. Members of NASA’s research staff also may be residents at ICASE 
for limited periods. 

The major categories of the current ICASE research program are: 

• Numerical methods, with particular emphasis on the development and analysis of basic 
numerical algorithms; 

• Control and parameter identification problems, with emphasis on effective numerical 
methods; 

• Computational problems in engineering and the physical sciences, particularly fluid 
dynamics, acoustics, and structural analysis; 

• Computer systems and software for parallel computers. 

ICASE reports are considered to be primarily preprints of manuscripts that have been 
submitted to appropriate research journals or that are to appear in conference proceedings. 
A list of these reports for the period October 1, 1989 through March 31, 1990 is given in the 
Reports and Abstracts section which follows a brief description of the research in progress. 


1 Presently, ICASE is operated at NASA Langley Research Center, Hampton, VA, under the National 
Aeronautics and Space Administration, NASA Contract No. NAS1-18605. In the past, support has been pro- 
vided by NASA Contract Nos. NASI- 18107, NAS1-17070, NAS1-17130, NAS1-15810, NAS1-16394, NAS1- 
14101, and NAS1-14472. 
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RESEARCH IN PROGRESS 


Saul Abarbanel 

The work on “spurious” oscillations due to non-characteristic imposition of boun ary 
conditions, reported in ICASE Report No. 90-16, has now branched out to become an 
investigation on possible interaction of inflow and outflow boundary conditions to produce 
oscillatory results with more than a single “spurious” frequency. This work, together with 
D. Gottlieb, is in progress and some of the new results wUl be reported on in the Third 
International Conference on Hyperbolic Problems to be held June 11-15, 1990 in Uppsa a, 

Sweden. , , , 

A new approach to accelerating convergence to steady state is also being stu le oge er 

with D. Gottlieb. The idea is to generalise the concept of local-time stepping by minimising 

the resulting residual (in a suitable norm that accounts also for boundary conditions). One 

then gets, even for linear problems, a distribution of Courant numbers that depends on the 

solution at the previous “time-step”. The method was tried on a model one-dimensional 

Burgers problem (with a stiff source term and a highly stretched mesh). The convergence 

rate in terms of the number of iterations required to reduce the residual by 5 orders of 

magnitude, showed an order of magnitude improvement (typically 4 vs. 55 iterations) m 

comparison to standard local-time stepping methods wherein the local Courant number is 

determined a priori from linear stability theory. This work will also be presented at the 

Uppsala conference. 


H. T. Banks and Fumio Kojima 

We are continuing our investigations (in collaboration with W. Winfree of the Instru- 
ment Research Division, LaRC) on inverse problems arising in thermal testing of materia s 
in space transportation systems and passenger airplanes. A spline-based parameter es i 
mation technique has been developed using techniques and ideas related to the ‘met o o 
mapping’. The ideas developed were successfully tested with the exper, mental data from 
IRD laboratories at NASA Langley. We are currently pursuing investigat.ons to further 
refine our methods with regard to ease and accuracy in detection and characterisation of 
flaws with thermal based methods. Furthermore, our results are now extended to the case of 
more complicated structures such as might occur in airplanes. The domain decomposition 
technique as well as the ‘method of mapping' plays an important role in our approach. 
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H. Scott Berryman, Charles Koelbel, Piyush Mehrotra, and Joel Saltz 

We are developing techniques which aUow a nonshared memory multiprocessor (such as 
ICASE’s iPSC/860) to be programmed at a higher level than is currently possible. Our 
research explores the efficient implementation, on distributed machines, of programs involv- 
ing accesses to globally defined data structures. It is natural to use operations on globally 
defined data structures to express data parallel algorithms typically used in scientific appli- 
cations. We have placed particular emphasis on examining optimizations that are suited to 
parallel loops over distributed globally defined arrays. 

Our work can be broken down into four interrelated areas: handling regular problems, 
handling irregular problems, run-time support, and compiler implementation. 

In all cases, our approach is to transform parallel do loops over globally defined data 
structures into code that executes on each processor. The code actually executed by each 
processor in a distributed memory machine can directly reference elements of global data 
structures that are stored on that processor. Data elements that are stored on other proces- 
sors must first be communicated by message passing. For efficiency, these nonlocal elements 
must be prefetched before the actual computation and stored in local memory. Two tech- 
niques are needed to deal with this prefetching: one for regular problems and one for irregular 
problems. 

In regular problems, such as solving PDEs on regular meshes and dense linear algebra 
applications, the pattern of nonlocal data accesses can be specified analytically at compile 
time, based on analysis of the program. In this case, it is possible to produce very efficient 
code by using this specification to generate the appropriate send and receive statements. We 
are pursuing this aspect of the work by deriving the necessary analytical formulas and by 
incorporating them into a compiler. 

In irregular problems, such as solving PDEs on unstructured meshes and sparse matrix 
algorithms, the communications pattern depends on the input data. This typically arises 
due to some level of indirection in the code. In this case, it is not possible to predict 
what data must be prefetched. We deal with this lack of information by transforming the 
original parallel loop into two constructs which we call the inspector and executor. During 
program execution, the inspector examines the data references made by a processor and 
calculates what off-processor data needs to be fetched and where that data will be stored 
once it is received. The executor loop then uses the information from the inspector to 
implement the actual computation. We are pursuing this aspect of the work by supplying user 
subroutines to implement inspector/executor pairs in existing languages and by generating 
the inspector/executor code in a compiler. We have described some of the facets of 
this work in “A Scheme for Supporting Automatic Data Migration on Multicomputers” 
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by Mirchandaney, Saltz, Mehrotra, and Berryman and in “Parallel Loops on Distributed 
Machines” by Koelbel, Mehrotra, Saltz, and Berryman. Both of these papers will appear 
in “Proceedings of the Fifth Distributed Memory Conference,” Charleston, South Carolina, 
April 1990. 

For run-time support we are developing a set of primitives (PARTI - Parallel Automated 
Runtime Toolkit at ICASE) which will be callable from FORTRAN or other languages. 
The PARTI system is a set of primitives that can be used to implement a wide range of 
scientific algorithms on distributed memory machines. These primitives support various run- 
time operations required by programs that need to make use of an embedded shared name 
space on a distributed machine. In addition, the primitives are used to study compile-time 
techniques that can be incorporated into a compiler for distributed memory architectures. 

The PARTI primitives are initialized by specifying a mapping into distributed memory 
for each globally defined multidimensional array. The primitives include procedures that 
allow one to scatter and gather array elements to the distributed memory. Primitives also 
are implemented that support several different methods of caching and retrieving copies of 
off-processor data. 

The PARTI system is organized into two levels. The lower level supports memory op- 
erations such as scatters and gathers across processors. The higher level binds mapping 
information to distributed arrays and uses this information to call the lower level primi- 
tives. The primitives allow the storage of information about memory access patterns so that 
memory operations with the same address calculations need not repeat these calculations. 
Send/receive schedules are generated for memory operations. The schedules may be stored 
and reused as well. This is particularly important in iterative codes. A preliminary descrip- 
tion of these primitives can be found in “Execution Time Support for Scientific Programs 
on Distributed Memory Machines” by Berryman, Saltz, and Scroggs; this will appear in 
the Proceedings of the SIAM Conference on Parallel Processing for Scientific Computing, 
Chicago, Illinois, December 1989. 

In the compiler implementation area, we have designed the Kali language as a testbed 
for implementing these techniques. The analysis for both regular and irregular problems is 
incorporated into the compiler. We are currently testing and enhancing the compiler and 
expect to eventually use it for realistic problems. 


Shahid Bokhari 

The multigrid method, an efficient procedure for solving partial differential equations, has 
very poor utilization when implemented on a massively parallel processor. We are exploring 
techniques for pipelining several parallel multigrid algorithms. By carefully organizing our 
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computations, we are able to achieve very significant increases in utilization. An ICASE 
report describing this work will be completed soon. 

A second effort is focused on the development of a network flow model for load balanc- 
ing in circuit-switched multicomputers. In multicomputers that utilize circuit switching or 
wormhole routing, communication overhead depends largely on link contention - the vari- 
ation due to distance between nodes is negligible. This has a major impact on the load 
balancing problem. In this case, there are some nodes with excess load (sources) and others 
with a deficit of load (sinks), and it is required to find which sources should transmit to which 
sinks so as to avoid contention. The problem is made complex by the fact that routing on 
currently available machines is hardwired. (On the Intel iPSC-2 and iPSC/860, the ‘e-cube’ 
algorithm is used; the Symult 2100 employs a ‘row-column’ strategy.) Thus the user can 
control only which nodes communicate but not how the messages are routed. 

Network flow models of message flow in the mesh and the hypercube have been developed 
to solve this problem. The crucial property of these models is the correspondence between 
minimum cost flows and correctly routed messages. To solve a given load balancing problem, 
a mi nimum cost flow algorithm is applied to the network. This permits us to determine 
efficiently a maximum contention free matching of sources to sinks which, in turn, tells us 
how much of the given imbalance can be eliminated without contention. An ICASE report 
describing this research is nearing completion. 

Experiments have been conducted on the Intel IPSC/860 hypercube in order to evaluate 
the overhead of interprocessor communication. It is demonstrated that (1) contrary to 
popular belief, the distance between two communicating processors has a significant impact 
on communication time, (2) edge contention can increase communication time be a factor of 
more them 7, and (3) node contention has no measurable impact. 


John Burns 

We are completing an ICASE report on the application of linear quadratic regulator 
theory to the problem of stability enhancement for the Burgers’ equation. In particular, 
optimal feedback gain operators were computed by minimizing an energy functional based 
on a linearization of the nonlinear model. These linear feedback control laws were then 
applied to the nonlinear Burgers’ equation and several theoretical and numerical results 
were obtained for the nonlinear closed-loop system. A computational scheme was developed 
to compute the optimal functional gains and several numerical experiments were preformed. 
The following important observations were made: 

1. The functional gains strongly depend on the location of the “sensors” used in the cost 
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functional. In the case where the cost functional contained “averaged gradients”, the 
functional gains were discontinuous at those points where the gradients were measured. 

2. The closed-loop system does decay at a faster rate than the open-loop (Burger’s equa- 
tion) and for large Reynolds number the steep gradients were smoothed by full linear 
state feedback. 

3. The LQR gains computed on low Reynolds number linearized models also enhanced 
stability for the nonlinear problem over a wide range of Reynolds numbers. Although 
the large gradients were still smoothed, performance was decreased. However, the 
closed-loop system was more stable than the open-loop. 

These results provide some insight into the more complex problem of fluid flow control. 

We axe continuing our efforts on active control of several fluid/structure problems. 


Richard Carter 

The interaction of structural design with control law formulation is an important topic 
in the construction of large scale, flexible space structures. By formulating the design pro- 
cess as a numerical optimization problem using engineering design goals as merit functions 
and constraints, improved structures and control laws can be determined simultaneously. 
The optimization problems arising from this approach are typically very computationally 
intensive, and are often highly nonlinear, nonconvex, and poorly conditioned, so that robust 
and efficient solution algorithms are imperative. In conjunction with the Control-Structures 
Interaction (CSI) group at LaRC, simulated model structures are being examined to develop 
methodology and to investigate features of optimization problems of this type. A variety 
of solution techniques are being applied to an Earth Pointing System (EPS) simulation 
model in order to identify the most promising approaches. Early difficulties associated with 
jump discontinuities and discontinuous gradients in some proposed design constraints/ goals 
have been alleviated by substituting carefully selected, smooth merit functions provided 
by Peiman Maghami of the CSI-ADM team. Several new approaches for self-scaling the 
perturbation level in finite difference gradient calculations are also being pursued. 

Another topic under investigation involves the robustness of optimization algorithms 
when gradient computations are inexact. Trust region methods and linesearch methods are 
the two major approaches to creating robust algorithms. Previous theoretical work and 
extensive numerical testing have shown algorithms based on the trust region methodology to 
be exceptionally robust with respect to gradient errors, while simple examples exist to suggest 
that linesearch algorithms are potentially very fragile with respect to even small errors. 
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Nevertheless, a recent software study at Cornell found two implementations of the linesearch 
approach that were competitive when finite difference gradients were used. I am currently 
examining these results to determine whether they are an artifact of the experimental design 
of the study (which compared different finite difference routines rather than specifically 
addressing trust region versus linesearch methodologies), or whether they reveal unknown 
principles about the behavior of linesearch methods. 


Leon M. Clancy 

ICASE has begun the process of converting its Sun network to SPARC based workstations 
running the 4.1 version of the operating system. In addition, two Silicon Graphics Personal 
Iris workstations have been integrated into the network. Other equipment upgrades include 
the installation of additional memory for SPARCstations and the file servers. 


Thomas W. Crockett 

The iPSC/2 hypercube computer was upgraded to an iPSC/860, with additional I/O 
capacity added. The system was checked out and made available to users. Accounting 
procedures were established to monitor processor and disk utilization. Benchmarks show 
that the upgraded system ranges from 5 to 35 times faster for double- precision floating- 
point computations. Substantial improvements are expected as the compilers mature. New 
releases of the X Window System and Andrew were installed on the ICASE Suns. A set of 
default configuration files were built to make it easier for novice users to get started with X 
and Andrew. 


Naomi Decker 

Research is proceeding along two paths with a common goal: an efficient and robust 
parallel multigrid algorithm for solving transonic and hypersonic compressible flow problems. 
The first path is the development of an good sequential multigrid code using semicoarsened 
multigrid algorithms with implicit relaxations. Past research in robust multigrid algorithms 
with J. VanRosendale showed that the coarsening of the grid in only one direction, and using 
an implicit method in the orthogonal direction(s), is an effective and efficient strategy for 
solving elliptic partial differential equations. The insensitivity of these methods to grids with 
high aspect ratio and to directionally-biased data dependencies has lead us to extend these 
methods to solve the non-elliptic systems arising in computational fluid dynamics. We wish 
to compare the efficiency of these new multigrid methods with easily parallelizable, explicit 
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Runge-Kutta, fully coarsened, multigrid algorithms. The Runge-Kutta multigrid algorithms, 
developed by A. Jameson for transonic calculations, can also be used for hypersonic problems, 

as shown in work with E. Turkel. # 

For the Euler or Navier-Stokes equations, the use of implicit time stepping requires the 

inversion of block tridiagonal systems, which can be efficiently solved using the block version 
of the Thomas algorithm for Gaussian elimination. However, since the implicit methods 
involve global data dependencies, their efficient parallel implementation is considerably more 
difficult and non-intuitive than for the explicit methods. Despite the extra effort involved, 
the potential for increased robustness justifies a careful examination of the implicit methods 
for parallel architectures. Thus, the other path is the parallelization of implicit schemes. We 
have considered various partitioning and scheduling schemes for alleviating the effects of t e 
global data dependencies. The communication and the computation aspects of these methods 
have been analyzed in the overall context of the solution process, including, for example, the 
effects of the various other data dependencies resulting from the boundary condition and 
residual computations. Performance of the proposed methods are being verified on .the 
Victor multiprocessor system, a message passing architecture developed at the IBM, 
Watson Research Center. A report summarizing these results is in preparation. 


Peter Duck , 

Work following on from ICASE Report Nos. 89-19 and 90-14 continue, on the temporal 

stability of axisymmetric supersonic boundary layers. Most recently the effect of wall sur ace 

cooling (and heating) has been investigated. This indicates that, just as ,n the planar case 

(which has been the subject of much investigation over the years), cooling genera y (i) 

stabilizes the so-called “first mode” of instability, but (ii) de-stabilizes the “second mo e 

of instability. Curvature is found to generally play a stabilizing role on the inv.scid stability 

of the Sow. In connection with this work, a number of computations are being carried out 

jointly with Dr. M. G. Macaraeg (Fluid Mechanics Division, LaRC) involving the full (sma 

disturbance) equations, with a view to comparing with a number of asymptotic theories, 

including those in ICASE Report Nos. 89-19 and 90-14. 

Work is also continuing with R. Bodonyi on 2D receptivity problems, in particular inves- 
tigating the interaction between a suction slot and a small amplitude unsteady freestream 
disturbance. These computations lead to a prediction of the amplitude of resulting Tollmien- 

Schlichting waves. . , 

Very recently, this type of work has been extended to three- dimensions. In particular the 

effect of a three-dimensional surface distortion inside a boundary layer on a small amplitude, 

unsteady, freestream disturbance is currently being computed. 
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Thomas Eidson 

In collaboration with M. Y. Hussaini and T. A. Zang (Fluid Mechanics Division, LaRC), 
the filtering analysis of a direct simulation (DNS) of the turbulent Rayleigh-Benard problem 
has been completed and a paper is currently being written. Briefly, a large eddy simulation 
(LES) type filter was used to split the velocity and temperature fields calculated in the DNS 
into a larger scale or filtered part and a smaller scale or fluctuation part. By varying the the 
filter width, information for a wide range of flow scales can be ascertained. The components 
of kinetic energy equations for the filtered and fluctuation velocity fields were computed and 
studied. One of the important results from this study was to quantify the magnitude of the 
thermal and “cascading” production of kinetic energy at different flow scales. Also the cal- 
culation of the diffusion, particularly for the small velocity scales, provided new information 
about turbulent flows. In addition, two LES models were evaluated by determining how well 
they predicted “cascading” production as compared to the direct calculation of the produc- 
tion from the simulation results. The scale similarity model gave much better predictions 
than the Smagorinsky model. 

A three-dimensional, incompressible, homogeneous shear flow code developed by T. A. 
Zang (Fluid Mechanics Division, LaRC) was studied for possible parallel optimization. This 
is a typical pseudo-spectral flow simulation code used at NASA and ICASE. Several parallel 
versions of the code were developed: 

1. small granularity, micro-tasking, 

2. large granularity, micro-tasking, 

3. large granularity, macro-tasking, 

4. large granularity, macro- tasking, minimum tasks. 

Granularity is defined as the amount of code that is executed in parallel at one time (before 
returning to the sequential controlling code). The small granularity version was over 10% 
faster than the other versions. This is possible due to the fact that the large granularity codes 
can have more code that is executed (redundantly) in each parallel process. The average 
speedup was 3.7 for the small granularity code and averaged 3.2 for the other versions. The 
study was done on a 4 processor CRAY-2. 


Gordon Erlebacher 

Work begun on compressible turbulence in collaboration with S. Sarkar, H. Kreiss, and M. 
Hussaini has been completed and reported in ICASE Reports 89-79 and 90-15. The objective 
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of this work was to understand the different types of turbulent structures as a function of 
initial conditions in compressible turbulent flows when the fluctuating Mach number is small. 
Several regimes were identified and explained for homogeneous isotropic turbulence when 
the fluctuating Mach number is small: a quasi incompressible regime, an acoustic regime, 
dominated by acoustic waves, and a regime where shocks are present. Studies are currently 

underway to study homogeneous compressible shear flow. 

Expertise was developed in the use of the new symbolic manipulation language, Math- 
ematica. Mathematica was then applied to compute the characteristics of several stability 
problems in incompressible and compressible flows. These problems include primary insta- 
bility, secondary instability, and weakly non-linear stability theory. The work related to 
weakly non-linear wave interaction in incompressible flows was done in collaboration with 
Bart Singer (High Technology Corporation, LaRC). New symbolic algorithms were developed 
which generate all the matrix coefficients of the discrete problem, along with over 8000 lines 
of fortran code in under 15 minutes on a personal Iris workstation. This execution time is at 
least two orders of magnitude faster than the brute force approach previously adopted with 
Macsyma. These symbolic manipulation programs, along with sophisticated linear stability 
programs allow users to quickly perform stability experiments using different sets of partial 
differential equations and different coordinate systems in a matter of hours instead of days. 

Side by side with the previous item, Herbert’s linear stability code linear.x was modified 
to provide a nice graphic interface. Linear.x allows users to quickly test out the stability 
properties of various problems with very little effort. Mathematica code has been written 
to interface directly with linear.x. The user inputs into Mathematica a set of differential 
equations, and in a matter of minutes, matrix coefficients have been computed and properly 
included into the code. To complement this effort, I have written a visual interface which 
allows linear.x to run on a Cray supercomputer, and display spectra and eigenfunctions 
on a local Iris workstation. The graphic interface permits the user to interactively pick 
eigenvalues off the screen, which are then sent back to the Cray. Then, in quasi real-time, 
the eigenfunction is computed via local procedure and returned to the Iris to be plotted. In 
this way, the user is freed from the burden of manipulating data files to ultimately visualize 

results. 


David Gottlieb 

Together with W. S. Don (Brown University) Mid J. Crawford (Brown University), we 
continue our research on flows past circular cylinders. We use the 2D code to study some 
fundamental features of the flow such as the critical Reynolds number and the effect of the 
actual body shape on the behavior of the flow. A 3D code is in its completion stages. 
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Together with S. Abarbanel we are studying the role of boundary conditions in creating 
numerical spurious frequencies. This problem has its origin in the study of flow past bluff 
bodies. We have shown that both the inflow and the outflow boundary conditions can create 

this type of frequency and we are studying the interaction between the inflow and outflow 
boundary conditions. 

We are developing essentially non-oscillatory spectral schemes for shockwaves calcula- 
tions. We have developed one side filters, e.g., filters that yield spectral accuracy in one side 
of the shock even at the shock itself. These filters are being tested and applied to two flows. 
We have shown that these methods can be applied successfully to interactions of shocks and 
density waves. We are trying to simulate shock- vortex interactions (with W. S. Don and A. 
Solomonoff (Brown University)). We continue our study in parallel spectral methods. The 
mam idea here is to use domain decomposition techniques and to use known conserved quan- 
tities in order to completely decouple the sub-domain. The runs are done on a hypercube 
architecture. 

Together with C. L. Streett (Fluid Mechanics Division, LaRC), we have developed accu- 
rate quadrature formula to substitute for the classical Clanshaw-Curtiss formula. We have 

shown that our formula is more robust when applied to spectral simulations of incompressible 
flows. 


M. YousufF Hussaini 

The research program deals with instability, transition and turbulence in simple canonical 
flows. Specifically, the linear stability of supersonic/hypersonic wedge and conical flows is 
being studied taking into account the presence of shocks. This study involves the assumption 
of parallel flow. It is planned in the next step to use parabolized stability equations which 
will include the nonparallel effects. The linear stability of compressible attachment-line 
boundary layer is also under study. The point of interest is the comparison of results from 
the two different approaches - the classical linear parallel flow stability theory and the triple- 
deck theory. In another study, the evolution of the compressible Couette flow from laminar 
to fully developed turbulent state is being investigated. The focus is the direct numerical 
simulation. Even so some aspects of the problem amenable to mathematical analysis are 
being studied. This research program is carried out in collaboration with G. Erlebacher and 
M. R. Malik (High Technology Corporation, LaRC). 


Thomas L. Jackson and Chester E. Grosch 

Our work focuses on the hydrodynamic stability of compressible high speed reacting free 
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shear flows, directly applicable to the study of supersonic diffusion flames in the context of 
scramjet engines. Our current work in this area includes the study of both near and far 
wakes, three-dimensional mean flows, and more realistic kinetics using quasi-global models 
for the chemistry. In tackling problems in this area, we have employed a combination of 
asymptotics and numerics to reduce the complex problems to model problems, thus isolating 
key physical effects for analysis. 


Peter A. Jacobs 

In recent years, shock-tunnels and expansion-tubes have been used to obtain experimental 
data at hypersonic flight speeds (3 < v < 7km/ s) for the National Aero-Space Plane project. 
The chemistry of the air is important at these speeds and much of the physical and numerical 
modelling used in Computational Fluid Dynamic (CFD) simulations has not been fully 
validated. We have been involved with the experimental side of this work, and now we are 
becoming involved with CFD simulations of the pulse facilities themselves and the generic 
scramjet models that have been tested in these facilities. 

We have also been studying the powder flow problems in the manufacturing process of 
Dynamic Powder Consolidation (DPC). In this process, a metal powder is compacted into 
a solid component by the passage of a strong shock wave. DPC has the ability to achieve 
compacts with close to theoretical densities and can do so without affecting the microstruc- 
ture of the bulk of the metal grains. We are modelling the powder as a compressible fluid 
with a nonstandard equation of state, and we are using a second-order Godunov method to 
compute the powder flow. 


Charles R. Johnson 

Among primary research projects completed during the period were a complete sparsity 
analysis of the QR factorization with D. Olesky and P. van den Driessche (University of 
Victoria, British Columbia) and a refined Gersgorin-type lower bound for the smallest sin- 
gular value with A. Hoffman(IBM Watson Research Center). In the former, given only the 
sparsity pattern of a matrix with independent columns, the zero patterns of both factors 
Q and R are predicted with algorithms that can be cheaply implemented. The results are 
exact in the sense that every position predicted to be nonzero will be nonzero for some (in 
fact, almost every) matrix with the original pattern, and every zero is necessary. This refines 
and completes work by several authors interested in algorithm design and storage allocation. 
Necessary zeros in Q or R can occur in a variety of subtle ways. 

In addition, a lengthy survey of matrix completion problems was completed for a volume 
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published by the IMA at Minnesota. 


Fumio Kojima 

Work is continuing on the development of parameter estimation techniques based on the 
boundary integral equation approach. Domain identification problems related to an electrical 
impedance tomography has been studied. Those are formulated as a parameter estimation 
problem for an integral equation of the second kind. We developed a spline based technique 
using the collocation method. Numerical experiments indicate the method is attractive and 
convergence arguments for the proposed algorithm are currently under study. 


Jacques Liandrat and Thomas W. Eidson 

Work is continuing on the definition and implementation of numerical algorithms for the 
resolution of partial differential equations based on a spatial wavelet approximation. Nu- 
merical results have been obtained in the case of the ID regularized Burgers equation with 
periodic boundary conditions. They show that algorithms based on the multiscale and lo- 
cal approach provided by the wavelet decomposition can be very efficient when a solution 
of a partial differential equation exhibits local behavior (rapid variations for instance) that 
requires locally very fine resolution. Different fields are now under investigation: the numer- 
ical analysis of such methods and their generalisation to non-periodic and multi- dimensional 
problems. Close to this numerical program is the use of wavelets as a tool to extract sub- 
stantial physical information from turbulence data. The study of the results of different 
numerical direct simulations of turbulence or transition is currently underway. The com- 
bination of these two approaches should lead to numerical simulations of turbulence on a 
wavelet basis. 


Dimitri Mavriplis 

Work is continuing on the use of unstructured triangular meshes for solving turbulent 
viscous flow problems about arbitrary configurations in two dimensions. In previous work , a 
method for generating highly stretched triangular meshes for accurate shear-layer resolution 
(ICASE Report No. 88-47) has been developed, and a finite-element Navier-Stokes solver 
has been devised and validated (ICASE Report No. 89-11). Hence, the present work has 
concentrated on refining the basic mesh generation and flow solution algorithms, as well as 
on the implementation of a suitable turbulence model, and the incorporation of adaptive 
meshing techniques in order to enhance the accuracy and efficiency of the method. 
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An algebraic turbulence model has been developed for use on unstructured meshes. Based 
on the Baldwin-Lomax model, this algebraic model makes use of background structured mesh 
lines. A local structured mesh is constructed in the vicinity of each solid wall boundary and 
wake line. At each time step within the solution of the flow on the global unstructured 
mesh, the current flow variables are interpolated onto these background structured meshes. 
The turbulence model is then executed on these meshes, and the resulting eddy viscosities 
are interpolated back onto the global unstructured mesh where they interact with the flow 
field solution. When adaptive meshing is invoked, the global unstructured mesh is refined 
by adding new points in regions of high flow gradients and performing local retriangulation. 
The turbulence background meshes are also refined by adding new normal mesh lines corre- 
sponding to each newly inserted unstructured mesh boundary point, and by increasing the 
resolution of the existing mesh lines in regions of high flow gradients. After each mesh adap- 
tation stage, the coefficients for interpolation between the new unstructured mesh and the 
refined background turbulence meshes must be recomputed prior to resuming time-stepping 
on the flow-field equations. This method has been employed to compute the turbulent flow 
over a 4-element airfoil configuration. Good comparison between computed and experimen- 
tal wind tunnel (surface pressure) measurements has been observed. Such problems are of 
particular importance to the aircraft industry, with regards to the design and development 
of high-lift airfoil-flap systems for take-off and landing conditions. 

Further work is aimed at investigating different turbulence models and extending this 

methodology to three dimensions. 


Kirsten Morris 

Robustness of controllers designed using Galerkin approximations has been shown for 
bounded semigroup control systems. The extension to “Pritchard-Salamon” systems is 
straightforward. However, accelerometers, which are becoming increasingly important in 
the study of flexible structures, do not fall within this framework. Under certain conditions 
of regularity, we have shown that control systems with accelerometers are “well-defined” in 
a certain sense, and that an equivalence exists between frequency domain and time domain 
definitions of stability. This last result is critical for modern controller design methods, such 
as Hoc techniques, which depend upon frequency domain criteria. This work is continu- 
ing with H. T. Banks. The eventual goal is to establish conditions under which the use of 
Galerkin techniques is justified in control of systems with these sensors. 

Research is being conducted with J. Juang (Spacecraft Dynamics Division, LaRC) into 
control of large scale space structures using methods. Second-order systems with col- 
located actuators and sensors can be stabilized by any controller with a frequency response 
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in the right-half plane. We are trying to establish a design technique for more general con- 
figurations which is robust with respect to damping errors. A preliminary design for active 
control of acoustic plane waves has been obtained. This design will be implemented by R. 
Silcox (Acoustics Division, LaRC). 


David M. Nicol 

While the parallelization of discrete-event simulations has been a topic of interest for at 
least ten years, few formal performance analyses of such systems have been developed. We 
have developed stochastic models of parallel discrete event simulations that help to answer 
a number of fundamental questions: what is the optimal performance one can achieve?, how 
well can so-called “conservative” methods perform?, and under what circumstances can a 
conservative method achieve better performance than an “optimistic” method? 

In another effort, we have been developing new algorithms for mapping highly irregular 
grids onto parallel machines and have been developing graphical tools to aid the user in 
viewing the partitioning of these grids. 


Yuh Roung Ou 

Research continued on control of fluid-structure interaction problems. The basic approach 
is to formulate the problem by coupling various control laws with theories of the Navier- 
Stokes equations. 

Our current effort is devoted to simulating the problem of controlling unsteady flow 
around a circular cylinder. The initial focus of our investigation is on the computational 
aspects of the cylinder undergoing a constant speed of rotation. The development of flow 
patterns and histories of lift and drag on the cylinder are computed. The good agreement 
between these computational results and the experimental results of Coutanceau and Menard 
confirm their respective validity. Subsequent work will remove the restriction of constant 
speed of rotation and extend our development to model large variations in rotation. The 
main objective of this study is to achieve the desired lift enhancement by varying the rotation 
rate. In particular, we want to incorporate control theory to explore the various theoretical 
aspects of our research. 

The ultimate goal of this project is to develop both computational and theoretical meth- 
ods for modeling the interaction of control and fluid dynamics. The tools developed here will 
be used to investigate fundamental questions regarding control of separated flows by using 
different boundary control mechanisms. 
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Merrell L. Patrick and Mark Jones 

Development of algorithms for solving the generalized eigenvalue problem is continuing. 
An eigensolver, LANZ, based on the Lanczos method has been implemented and tested 
on the Convex, Cray Y-MP, and Cray 2 computers. The performance of a uni-processor 
version of LANZ has been compared with that of the eigensolver module, EIG, in the NASA 
Langley structural analysis testbed system running on the Convex and Cray 2 systems. For 
all problems tested the performance of LANZ is substantially superior to that of EIG. 

Our approach required the development of a routine for solving indefinite linear systems of 
equations. Such systems arise because shifting is required to improve convergence rates of the 
Lanczos iterations. A shifting strategy was developed for improving the performance of the 
algorithm and tested for robustness. A routine based on the Bunch-Kaufman algorithm for 
solving indefinite systems was implemented and shown to perform better than a LINPACK 
algorithm based on LU factorization. Three new algorithms, based on Bunch-Kaufman 
factorization, that take advantage of high performance architectures have been developed. 

A multiprocessor version of LANZ has been implemented and tested using the Force 
language on the Cray Y-MP computer. Preliminary results suggest three important activi- 
ties for the future for obtaining a more efficient implementation, particularly for very large 
problems. They are: 

1. The factorization phase of the Lanczos process should be as efficient as possible. We 
intend to add the option of sparse factorization to LANZ which currently uses variable 
band factorization. 

2. Since backsolving will not parallelize well, we need to reduce the number of backsolves 
in the Lanczos process. This can be accomplished by reducing the number of Lanczos 
iterations which calls for intelligent dynamic shifting. 

3. Iterative methods require less storage and parallelize more easily than direct elimination 
methods. The option of iterative methods will be added to LANZ if they prove to offer 
advantages over elimination methods. 

We will also begin to make an assessment on the appropriateness of other parallel architec- 
tures for implementing LANZ. 

This work is being done in collaboration with the NASA Langley CSM group. 
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Merrell L. Patrick and Terrence W. Pratt 

Methods for porting existing FORTRAN based applications codes to new parallel archi- 
tectures are the focus of this study. The NASA LAURA code (Langley Aerothermodynamic 
Upwind Relaxation Algorithm), developed by Peter Gnoffo (Space Systems Division, LaRC), 
is being studied as a test case. We are looking at both spatial and temporal partitioning 
methods for parallel execution of this code on various architectures, including the Intel 
iPSC/860. An execution time profile of the code has been obtained and an analysis of data 
partitioning and data flow from subroutine to subroutine is being carried out. 


Ugo Piomelli 

Further studies on the applicability of the large-eddy simulation approach to transitional 
wall-bounded flows have been performed. In large-eddy simulations (LES) of the Navier- 
Stokes equations, the large, energy- carrying scales of motion are accurately captured while 
the effect of the small scales is modeled. Thus, LES provides an attractive alternative to 
both the direct simulation, in which all scales of motion are resolved, and which require very 
massive computational efforts, and to the solution of the Reynolds averaged equations of 
motion, which require modeling of all scales of motion. The models used to parameterize the 
effect of the subgrid scales (known as residual stress models or subgrid-scale stress models) 
are simpler and more universal than the turbulence models used for the Reynolds-averaged 
Navier-Stokes equations, since the subgrid scales tend to be more isotropic and homogeneous 
than the large ones, and less affected by the boundary conditions. This technique has been 
successfully applied to many turbulent flows, but only recently have efforts been made to 
study transition to turbulence using LES. While during the linear stages, very little subgrid 
scale activity occurs, and subgrid-scale models should contribute little to the Reynolds stress, 
during laminar breakdown, significant amounts of small scales are generated, and inclusion of 
a model greatly improves the accuracy with which statistical quantities of engineering interest 
can be predicted. Since standard subgrid-scale stress models damp the growth of the velocity 
perturbation at the early stages of transition, modifications were developed to obviate this 
difficulty. The results of the calculation of transition to turbulence in a plane channel indicate 
that such models allow one to predict the entire development of the perturbations from the 
linear stage through laminar breakdown into fully- developed turbulent flow. Current efforts 
are directed towards the development of more general models that better incorporate the 
flow physics. 
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Peter Protzel 

Work continued on the investigation of the fault- tolerance characteristics of artificial neu- 
ral networks (ANNs) and their applications in the area of Reliable Computing. Together 
with D. Palumbo (Information Systems Division, LaRC), we studied the application of an 
ANN to control the allocation of tasks to processors in a distributed multiprocessor sys- 
tem. The architecture of the multiprocessor system is based on the Software Implemented 
Fault-Tolerance (SIFT) Computer, which was developed by NASA as a highly reliable, recon- 
figurable computer system for flight control applications. In SIFT, the failure of a processor 
requires the isolation of this processor and a reconfiguration of the system with a reallocation 
of all tasks among the remaining processors. The reallocation has to observe hard real-time 
constraints and should be performed in a way that balances the load of the processors. A 
certain type of ANN that solves constraint optimization problems in real-time can be used to 
generate these new allocations. With the ANN being a critical component in this application, 
the fault-tolerance of the ANN itself becomes an important characteristic. By studying the 
performance of the ANN in the presence of simulated hardware faults, we could demonstrate 
the extreme fault-tolerance with a hardly discernible degradation even after multiple faults. 
These characteristics of high speed, fault-tolerance, and low weight and power consumption 
make a hardware implemented ANN interesting for this and similar applications in the area 

of real-time control. 

Another project in collaboration with C. Jeffries from Clemson University investigates 
a novel architecture of a high order ANN that functions as an auto-associative memory. 
The main feature of this high order model is the absence of spurious attractors and the 
guaranteed convergence to one of the stored memories. One application for this ANN is the 
decoding and error-correction of block codes transmitted over a noisy channel. Since the 
ANN converges to the closest attractor (code word) stored in memory given an arbitrary 
analog vector as its initial input, it is possible to perform soft-decision decoding on block 
codes with the corresponding increase in performance. In computer simulations, we studied 
the performance of the ANN in terms of post-decoding bit error rate versus signal-to-noise 
ratio for two exemplary block codes. The results in comparison with a conventional decoding 
algorithm show, for example, that the bit error rate for a (15,5) block code at 7dB signal-to- 
noise ratio can be decreased by two orders of magnitude. A hardware implementation could 
produce different integrated circuits for different block codes to be employed as real-time, 
soft-decision decoders in potentially many different communication systems. 


Joel Saltz and Ravi Mirchandaney 

Dependencies between loop iterations cannot always be characterized during program 
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compilation. Doacross loops typically make use of a priori knowledge of inter-iteration 
dependencies to carry out required synchronizations. We propose a type of doacross loop that 
allows us to partition iterations of a loop between processors without advance knowledge of 
inter-iteration dependencies. The method proposed for partitioning loop iterations requires 
us to carry out parallelizable stages of preprocessing and post-processing during program 
execution. We have performed experimental studies using this construct and will report our 
results in an upcoming ICASE Report. 


Joel Saltz and V. Venkatakrishnan 

A fully implicit scheme is used to solve the compressible Euler equations for transonic 
flows. An approximate Riemann solver is used to compute the numerical fluxes to second 
order accuracy in space. The large linear system arising at each time step from the lineariza- 
tion of a fully implicit scheme is solved by iterative means using a preconditioned conjugate 
gradient-like algorithm. Incomplete LU (ILU) factorization with no fill-in is used to pre- 
condition the system followed by the use of the Generalized Minimum Residual technique 
(GMRES). Two ways to achieve parallelism are tested: one which makes use of parallelism 
inherent in triangular solves and the other which employs domain decomposition techniques. 
The vectorization/parallelism in triangular solves is realized by the use of a reordering tech- 
nique (called wavefront ordering) wherein the triangular matrix is interpreted as a directed 
graph and the data dependencies are analyzed. The graph is analyzed once using the con- 
nectivity of the grid points and all the matrix operations are done in the context of 4 x 4 
blocks since we have four unknowns to a grid point. The factorization can also be done in 
parallel with the wavefront ordering. In the domain decomposition case, the wave front or- 
dering is used for getting vectorization within each sub-domain and, since the sub-domains 
are disjoint, the factorization and triangular solves axe done in parallel. The domain de- 
composition does introduce a sequential overhead. We compare the performances of two 
ways of partitioning the domain, viz. strips and slabs. Results on a Cray Y-MP are re- 
ported for an inviscid transonic test case. The performances of the linear algebra kernels are 
also reported. This work is reported in “Domain Decomposition in Aerodynamics” in the 
Proceedings of SIAM Conference on Parallel Processing for Scientific Computing, Chicago, 
Illinois, December 1989. 


Sutanu Sarkar 

We are engaged in the direct simulation and Reynolds stress modeling of turbulent flows. 
Most of this research addresses compressibility effects in high-speed turbulence. In collabo- 
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ration with G. Erlebacher and M. Y. Hussaini, we have performed simulations of isotropic 
turbulence at various Mach numbers and are presently engaged in the simulation of homo- 
geneous shear turbulence with a spectral collocation method. These simulations and an 
asymptotic theory (developed by H. 0. Kreiss in association with G. Erlebacher, M. Y. 
Hussaini and myself) are being used to assist in the generation of turbulence models. We 
have obtained an interesting result that compressibility leads to enhanced dissipation, and a 
model parameterizing this phenomenon has been developed and successfully applied to the 
supersonic shear layer. We are continuing the application of second-order closures to other 
high-speed flows in collaboration with L. Balakrishnan (Research Associate, Old Dominion 

University). 


Paul Saylor 

An ICASE report on iterative methods and initial value problems, written with R. Skeel 
of the University of Illinois, is nearing completion. Analysis in the paper suggests that if the 
linear systems in a stiff initial value problem are to be solved iteratively, then the Chebyshev 
iterative method would be appropriate. This method is implemented in the package Cheby- 
code of Ashby and Manteuffel. The differential algebraic system solver (DASSL) package of 
Petzold has been modified (at the University of Illinois and Lawrence Livermore National 
Lab) to call Chebycode. The resulting package will be applied to a discretization of the low 
frequency approximation of the transonic small disturbance equation. 

Another effort is focused on separating eigenvalues. The small interior eigenvalues for a 
generalized eigenvalue problem are difficult to compute. It is an unresolved question whether 
using a polynomial to separate the eigenvalues in conjunction with a Lanczos-style algorithm 
would be a good idea. If it is, then there exist algorithms and software, due to S. F. Ashby, 
for computing a polynomial of given degree that separates eigenvalues optimally. To explain, 
suppose that the interval [A b , A c ] is contained in a cluster of eigenvalues. Let fi £ [A fc , A c ] . 
Ashby’s algorithms, based on the Remez algorithm, compute a polynomial, w fc , of given 
degree, k, such that (i) u k is the minimax polynomial on the complement of [A fc ,A c ]; and 
(ii) tu fc (/i) = 1. The peak at fi separates the eigenvalues clustered around fi from all other 
eigenvalues. Some technical details are omitted; for example, the complement of [A t , A c ] 
means the complement within the spectrum. A careful definition of the complement makes 

oifc unique. 

Polynomial is the polynomial approximation of the delta function on the complement 
of [At, A c ]. Polynomial approximations to the delta function may also be computed using 
orthogonal polynomials in which case they are called kernel polynomials. I have previously 
discussed using kernel polynomials to separate eigenvalues. Aside from questions about the 
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superiority of the minimax norm, the advantages to the minimax polynomial approach is 
the existence of algorithms and their Fortran implementation. 

Ashby’s application was to the iterative solution of symmetric indefinite linear algebraic 
equations. In that case, the objective is to find a polynomial, u k , that is small on the positive 
and negative eigenvalues and takes on the value 1 at fi — 0. The particular value /j, — 0 is 
easily generalized in the algorithm. A description of the algorithms and the application is 
in the December, 1989, special issue of BIT on preconditioning. 


Paul Saylor and Jeffrey Scroggs 

This is a preliminary investigation into combining asymptotic analysis and waveform 
relaxation to form a numerical method with both high accuracy and parallel capability. 
Since waveform methods easily exploit medium-grain to small-grain parallelism, the focus 
will be on obtaining accuracy. 

In its application to CFD, the method of lines yields a system of ordinary differential 
equations by discretizing the spatial derivatives. The waveform technique is to decouple 
the ODEs, solve each equation, then employ functional iteration to obtain a solution to the 
system. Waveform relaxation has the advantages that different time steps can be used for 
different components, and that decoupling makes it trivial to distribute the system among 
multi-processors for a parallel solution. 

The role of asymptotic analysis is to adapt waveform relaxation to the special require- 
ments of CFD. As an example, asymptotic analysis may be used to formulate a conservative 
scheme in the presence of shocks. In addition, asymptotic analysis may be used to resolve the 
physically relevant and diverse scales that arise when simulating flow with smooth regions 
interspersed with boundary layers and shocks. 

The first phase of this project is to demonstrate the waveform technique on a conservative 
finite difference approximation to the fluid dynamic equations representing conservation of 
mass and conservation of momentum. We are implementing the technique using DASSL, the 
differential algebraic system solver of L. Petzold. This work is being done in collaboration 
with Fen- Lien Juang of the University of Illinois. 


Jeffrey Scroggs and Marc Garbey 

This research is the symbiosis of asymptotics and numerical analysis to generate domain 
decomposition methods that are both accurate and computationally efficient. 

We are developing a method that consists of multiple- stages. The first stage of the method 
is to obtain a first approximation by a standard conservative numerical method (e.g., the 
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Godunov method). Subsequent stages of the method consist of corrections. The correction 
differs depending on whether the problem is inviscid (hyperbolic) or viscous (parabolic). A 
shock-layer correction which minimize the effect of the viscosity is used for the inviscid case; 
whereas, the viscous profile of the shock is computed when the problem is viscous. The 
method is derived and justified via asymptotic analysis that treats the order of accuracy 
e = 0(At p , Ax q ) of the standard numerical method as a small parameter in an asymptotic 
expansion. 

A residual correction based on asymptotic analysis is utilized outside the zone of strong 
singularities. For inviscid problems, this residual correction can be understood as an antid- 
iffusion correction. For example, this correction is able to improve the Godunov scheme in 
the neighborhood of a weak singularity such as a jump of the first derivatives of the solution 
(corner-layers). 

The method has been implemented for a system of equations arising in fluid dynamic 
equations representing conservation of mass and conservation of momentum, and has demon- 
strated that capturing the shock-layer profile (for a slowly varying profile) can be performed 
in an efficient manner. 


Sharon S. Seddougui 

The spatial inviscid instability problem for Gortler vortices is being investigated jointly 
with Philip Hall. The aim of this work is to determine the fastest growing spatial mode for 
large Gortler numbers. For an incompressible fluid, it is found that these modes are localized 
in an asymptotically thin layer at the wall. The receptivity problem has also been addressed 
showing that 0(1) coupling coefficients between a surface perturbation and the induced 
vortex field are possible. The extension of this work to include the effects of compressibility 
is currently being studied. 


Chi- Wang Shu 

We are continuing our investigations of ENO (essentially non-oscillatory) finite difference 
schemes. One project is to obtain a user-oriented two-dimensional Navier-Stokes code which 
handles generalized coordinates via transformations and various boundary conditions and 
to use it to test problems including direct turbulence simulation and acoustic related phe- 
nomena. This project is coordinated with the Theoretical Flow Physics Branch of LaRC. 
Another project is to analyze the accuracy of ENO schemes and how to effectively use cen- 
tered schemes in the smooth part of the flow and use ENO in the rapid transition regions. 
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Charles G. Speziale 

Research has continued on the development of a compressible second-order closure model 
for high-speed turbulent flows in collaboration with S. Sarkar. This latest research has con- 
centrated on the development of models which are asymptotically consistent near walls and 
computationally robust. Some preliminary boundary layer calculations conducted with R. 
Abid (Vigyan, LaRC) have indicated that none of the existing second-order closure models 
have both of these desired properties. In order to properly account for anisotropies in the 
turbulent dissipation rate, a modeled tensor dissipation rate transport equation was recently 
developed in collaboration with T. B. Gatski (Fluid Mechanics Division, LaRC). A more 
simple algebraic model will be extracted from this transport model based on a local homoge- 
neous equilibrium hypothesis. Unlike the older models, this new algebraic model will account 
for nonlinear effects arising from the mean velocity gradients. The compressible dissipation 
model recently developed by Sarkar, Erlebacher, Hussaini, and Kreiss (ICASE Report No. 
89-79) will then be incorporated into this algebraic model to account for compressible effects. 

Research continues on the large-eddy simulation of transitional flows with U. Piomelli and 
T. A. Zang (Fluid Mechanics Division, LaRC). The most recent discussions have centered 
around the development of a one-equation subgrid-scale model to account for history effects. 
Results obtained previously were based on a rescaled Smagorinsky model which does not 
directly account for these effects which are known to be important in most transitional 
flows. New work on the large-eddy simulation of compressible turbulent flows - conducted 
in collaboration with G. Erlebacher, T. A. Zang, and M. Y. Hussaini — has also begun 
recently. 


Shlomo Ta’asan 

The development of efficient multigrid solvers for constraint optimization problems gov- 
erned by partial differential equations has continued. A new method for treatment of control 
problems governed by elliptic PDE’s has been developed. The case of finite dimensional con- 
trol has been considered so far. In that case, solution for the full optimization problem was 
reached with the same cost as solving the constraint equations about 2-3 times. Research 
is continuing in applying the techniques developed for identification problems in which the 
parameter space is finite dimensional (and not too large), e.g., impedance tomography. 

The same ideas axe being considered in the context of aerodynamics design problems 
where airfoils are to be calculated so as to meet certain design requirements. For example, 
to give pressure distribution in some flow conditions which are closest to a given pressure 
distribution. In the course of developing these ideas new Euler solvers are being developed 
which uses DGS relaxation schemes and perform equally well for all Mach numbers. 
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Another area of research is the development of new computational techniques in elastic- 
plastic problems. Here a set of non-linear equations which involve also weak dependence of 
the history of the loading are to be calculated. Our effort is focused on trying to perform 
most of the computation on coarse levels. Different continuation techniques (in the load 
parameters) are being tested to obtain ones which are best suited for multigrid algorithms. 

A new method for the calculation of several eigenfunctions of elliptic problems has been 
developed. Its computational complexity is O(qN) where q is the number of eigenfunctions 
involved and N is the number of unknowns for each eigenfunction. Previously developed 
algorithms had complexity of 0(q*q*N), because of the Ritz projection that is needed for 
such problems. Our new algorithm also uses a Ritz projection but it is done on the coarsest 
level so its cost is negligible, reducing the cost of the algorithm to a minimum. 


Eitan Tadmor 

We study the convergence rate of approximate solutions to the possibly discontinuous 
entropy solutions u{x,t) of scalar nonlinear conservation laws. In [ICASE Report No. 89-80] 
we have shown that by post-processing the small viscosity approximations u e , we can recover 
pointwise values of u and its derivatives with an error as close to e as desired. 

Our approach does not ‘follow the characteristics’ and, therefore, applies mutatis mutan- 
dis to other approximate solutions such as E-difference schemes. Together with H. Nessyahu 
we extend these results to E-stable approximations, Glimm-type schemes and we are working 
on the convergence rate of certain kinetic approximations discussed below. 

Together with B. Perthame (Universite de Mathematiques, France) we study the con- 
struction of multi-dimensional entropy solutions for nonlinear scalar conservation laws. To 
this end we employ kinetic solutions depending on the microscopic scale t > 0 which were 
constructed in ICASE Report No. 90-11. We show that the phase averages (as well as higher 
moments) of these kinetic solutions, form a compact family of approximate solutions which 
are strongly convergent to the ‘macroscopic’ quantities associated with the multi-dimensional 
conservation laws (e.g., fluxes, entropy, etc.), as the microscopic scale, e j 0. 

Moreover, we precisely identify the governing evolution model of the kinetic limit solution, 
as t j, 0. The governing model is shown to consist of linear transport equations with a phase 
space gradient of the sink terms which are responsible for the loss of entropy at the kinetic 
level. Using this we conclude that for multi-dimensional scalar conservation laws which 
satisfy a rather weak nonlinearity assumption, then the corresponding macroscopic solution 
operator is compact . 

Finally, together with D. Gottlieb, we study the stability of spectral approximations to 
scalar hyperbolic initial-boundary value problems with variable coefficients. Time is dis- 
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cretized by explicit multi-level or Runge Kutta methods of order < 3 (forward Euler time 
differencing is included), and we study spatial discretizations by spectral and pseudospectral 
approximations associated with the general family of Jacobi polynomials. We prove that 
these fully explicit TV- degree spectral approximations are stable provided their time-step, At, 
is restricted by the CFL-like condition, A t < Const. N~ 2 . We give two independent proofs 
of this result, depending on two different choices of appropriate L 2 -weighted norms. Our 
result confirms the commonly held belief that the above CFL stability restriction, which 
is extensively used in practical implementations, guarantees the stability (and hence the 
convergence) of fully- explicit spectral approximations in the non-periodic case. An ICASE 
report describing this research is nearing completion. 


Bram van Leer 

The work on the multi-dimensional Riemann Solver, in collaboration with C. Rumsey 
(Fluid Mechanics Division, LaRC), has led to a robust algorithm that treats shocks and 
shears in the proper frame and thereby considerably reduces their numerical spread. After 
implementation in a 2-D Euler airfoil code, the next step is incorporation of the new flux 
function in a Navier-Stokes code. The results of the Euler calculations have been summarized 
in an extended abstract submitted to the next AIAA Aerospace Sciences Meeting, January 
1991, Reno, Nevada. 


John Van Rosendale and Naomi Decker 

In recent work, we have been looking at the design of iterative methods for the linear 
systems arising in implicit spectral calculations. Our first target is the variable coefficient 
Helmholtz equation, a problem arising in the Ozawa formulation of the Navier-Stokes equa- 
tions and in stream- function/ vorticity ocean circulation models. 

Deville and Mund have explored the idea of using finite elements as preconditioners for 
spectral methods. They found that standard trilinear finite elements in three dimensions lead 
to a preconditioned spectral system having condition number 1.44, a major improvement over 
finite difference preconditioned spectral methods. We have improved this condition number 
to 1.25. The new discretization is formed by using linear finite elements in one dimension with 
a partially lumped mass matrix. One then tensors this discretization appropriately to yield 
the optimal 9-point discretization in two dimensions, and optimal 27-point discretization in 
three dimensions. The new discretization induces no extra work and retains the standard 
tensor product structure, allowing efficient treatment of variable coefficient problems without 
storage of multiple coefficient arrays. 
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We have also designed a Petrov-Galerkin discretization based on quintic splines, leading 
to a preconditioned system having condition number 1 .001, for smoothly varying coefficients 
To exploit this low condition number, we constructed a two level algorithm with the spectral 
method preconditioned by this spline discretization, and the latter, in turn, preconditioned 
by a 27-point finite element discretization. The resulting algorithm reduces the spectra 
residual by a factor of a million, through six solves of the finite element linear system, but 

requires only one spectral residual calculation. . , 

To solve the finite element systems occurring, we use, ZOOM, a fast semicoarsening ase 
multigrid algorithm. ZOOM has a spectral radius of .003 per v-cyde on the positive definite 
Helmholtz equation on a uniform grid. It’s convergence rate degenerates to about .05 on 
Chebyshev grids, but this is still fast enough that one v-cycle suffices for each solution o 

the finite element linear system above. 


John Van Rosendale and Piyush Mehrotra , . . 

In the last year, we have been designing a programming environment, called Kali, which 
provides a software layer supporting a global name space on distributed memory arch, ec 
tures In Kali, the computation is specified via a set of parallel loops, whic act on va ues 
in this global name space. The user provides annotations to specify distribution of arrays 
across the distributed memory, and to assign work to processors. In this way, the user 
control, aspects of the program critical to performance, such as data distribution and load 

balance, while relegating the details of message-passing to the compiler. 

We are currently studying the efficacy of this approach using a Fortran-like language, 
KF1 (Kali Fortran 1). KF1 contains most of standard Fortran, plus Kali-style construe s 
for expressing parallel algorithm, for distributed execution. While no KF1 compiler exists 
yet, we have been studying the efficacy of KF1 constructs using an existing compiler for 

Pascal-like version of Kali. . . 

Preliminary results suggest that KF1 constructs are quite effective for expressing regular 

algorithms, such as tensor product array computations, relying on static load balancing^ 
Work on several large numerical applications programs, including a multign a gon m an 
an ADI algorithm, suggests that programs written in KF! will obtain efficiency comparable to 
that of analogous programs expressed in lower level message-passing languages. We have also 
recently added constructs to the language for problems requiring dynamic load balancing, 
and have begun looking at irregular problems requiring run-time analysis. 
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step in an explicit time marching algorithm is 0{N~ 2 ) which, in many cases, is much below 
the time step dictated by the physics of the P.D.E. In this paper we introduce a new set of 
interpolating points such that the eigenvalues of the differentiation operator are 0(N) and 
the allowable time step is 0(N ). The properties of the new algorithm are similar to those 
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Laarhoven (1985) and point out the possible drawbacks of these algorithms. By including 

algorithm to 75 ) p ^ Jec f tlons in the variable metric updating we can generalize 7 Straeter’s 
algorithm to a family of parallel projected variable metric algorithms which do not suffer the 
above drawbacks and which retain quadratic termination. Finally we consider the numerical 

hAw°tV^ ai l Ce - 0f °f t k 16 ? be i f ° f tbe fanul y on several standard example problems and illustrate 
how the choice of the displacement vectors affects the performance of the algorithm. 


L ’’ A ’ Baylis ®> f n ? R- Krishnan: Numerical study of three-dimensional spatial in- 
stability of a supersonic flat plate boundary layer. ICASE Report No. 89-74. October fi 1Q8Q 
19 pages. Submitted to Physics of Fluids. ’ ’ 


is studied eh TbL°nhl 8P r Uy ^ owm g three-dimensionS waves in a supersonic boundary layer 
is studied The objectives are to validate our computer code by comparing with experiments 

and parallel and non-parallel linear stability theories for = 4.5, ^S/ft = 1550 and at 
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T t>00 = 100°F and then use the program to study nonlinear effects. The three-dimensional un- 
steady Navier-Stokes equations are solved by a finite difference method which is fourth-order 
and second-order accurate in the convection and viscous terms respectively, and second-order 
accurate in time. Spanwise periodicity is assumed. The inflow disturbance is composed of 
eigenfunctions from linear stability theory. Computed results for small amplitudes of this 
disturbance agree well with linear theory and experiment for several frequencies. By increas- 
ing the amplitude of the inflow disturbance, nonlinear effects in the form of a relaxation type 
oscillation of the time signal of pu are observed. 


Brewer, Dennis W., John A. Burns, and Eugene M. Cliff: Parameter identification for an 
abstract Cauchy problem by quasilinearization. ICASE Report No. 89-75, October 13, 1989, 
35 pages. Submitted to the Quarterly Journal of Applied Mathematics. 


A parameter identification problem is considered in the context of a linear abstract 
Cauchy problem with a parameter- depen dent evolution operator. Conditions are investigated 
under which the gradient of the state to a parameter possesses smoothness properties which 
lead to local convergence of an estimation algorithm based on quasilinearization. Numeri- 
cal results are presented concerning estimation of unknown parameters in delay-differential 
equations. 


Nicolaides, R. A.: Direct discretization of planar div- curl problems. ICASE Report No. 89-76 , 
October 23, 1989, 31 pages. Submitted to SIAM Journal. 


A control volume method is proposed for planar div-curl systems. The method is inde- 
pendent of potential and least squares formulations, and works directly with the div-curl 
system. The novelty of the technique lies in its use of a single local vector field componexit 
and two control volumes rather than the other way round. A discrete vector field theory 
comes quite naturally from this idea and is developed in the paper. Error estimates are 
proved for the method, and other ramifications investigated. 


Bassom, Andrew P. and Sharon 0. Seddougui: A note concerning the onset of three di- 
mensionality and time dependence in Gortler vortices. ICASE Report No. 89-77, Octo- 
ber 31, 1989, 28 pages. Submitted to Journal of Fluid Mechanics. 


Recently Hall & Seddougui (1989) considered the secondary instability of large amplitude 
Gortler vortices in a growing boundary layer into a three-dimensional flow with wavy vortex 
boundaries. They obtained a pair of coupled, linear ordinary differential equations for this 
instability which constituted an eigenproblem for the wavelength and frequency of this wavy 
mode. Investigations into the nonlinear version of this problem by Seddougui & Bassom 
have revealed several omissions in the numerical work of Hall & Seddougui; these issues 
are addressed in this note. In particular, we find that many neutrally stable modes are 
possible; we derive the properties of such modes in a high wavenumber limit and show that 
the combination of the results of Hall & Seddougui and our modifications lead to conclusions 
which are consistent with the available experimental observations. 
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Zang, Thomas A. and M. Yousuff Hussaini: Multiple paths to subharmonic laminar break- 
down in a boundary layer. ICASE Report No. 89-78, October 31, 1989, 14 pages. Physics 

Review Letters, in press. 


Numerical simulations demonstrate that laminar breakdown in. a boundary layer in- 
duced by the secondary instability of two-dimensional Tollmien-Schlichting waves to three- 
dimensional subharmonic disturbances need not take the conventional lambda vortex/high- 
shear layer path. 


Sarkar, S., G. Erlebacher, M. Y. Hussaini, and H. 0. Kreiss: The analysis and modeling of 
dilatational terms in compressible turbulence. ICASE Report No. 89-79 , December 6, 1989, 
33 pages. Submitted to Journal of Fluid Mechanics. 


It is shown that the dilatational terms that need to be modeled in compressible turbulence 
include not only the pressure-dilatation term - the compressible dissipation. The nature of 
these dilatational terms in homogeneous turbulence is explored by asymptotic analysis of 
the compressible Navier-Stokes equations. A non-dimensional parameter which characterizes 
some compressible effects in moderate Mach number, homogeneous turbulence is identified. 
Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, 
and their results are found to be in agreement with the theoretical analysis. A model for 
the compressible dissipation is proposed; the model is based on the asymptotic analysis 
and the direct numerical simulations. This model is calibrated with reference to the DNS 
results regarding the influence of compressibility on the decay rate of isotropic turbulence. 
An application of the proposed model to the compressible mixing layer has shown that the 
model is able to predict the dramatically reduced growth rate of the compressible layer. 


Tadmor Eitan: Local error estimates for discontinuous solutions of nonlinear hyperbolic 
equations. ICASE Report No. 89-80 , December 6, 1989, 21 pages. Submitted to Communi- 
cations on Pure & Applied Mathematics. 


Let u(x,t) be the possibly discontinuous entropy solution of a nonlinear scalar conser- 
vation law with smooth initial data. Suppose u e (x, t) is the solution of an approximate 
viscosity regularization, where t > 0 is the small viscosity amplitude. We show that by 
post-processing the small viscosity approximation u e , we can recover pointwise values of u 

and its derivatives with an error as close to e as desired. , 

The analysis relies on the adjoint problem of the forward error equation, which m this 
case amounts to a backward linear transport equation with discontinuous coefficients. The 
novelty of our approach is to use a (generalized) E-condition of the forward problem in order 
to deduce a IF 1 ’ 00 energy estimate for the discontinuous backward transport equation j this, 
in turn, leads us to e-uniform estimate on moments of the error u e u. 

Our approach does not 'follow the characteristics’ and, therefore, applies mutatis mutan- 
dis to other approximate solutions such as E-difference schemes. 
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Hall Philip- Gortler vortices in 3 romn, boundary layers: The leading edge ereceptiv, proi- 
“m, linear growth, and the nonlinear breakdown stage. ICASE Report No. 89J1 , Decern 

bei^H, 1989, 74 pages. To appear in Mathematica. 


Gortler vortices are thought to be the 

importance. In this paper a govern this g growth and parallel flow theories do 

“g^ the' paraUel C'Zries h“ e some 

considered. It. is found that the mo. t e^d^tth^e “Ible experimental 
neutral curve is det J™ n ^f ffe ^\\ C e g imes growt h of nonlinear vortices is also given. Again 
data. A discussion of the different regimes 8 nonoarallel effects are dominant, 

it is shown that, unless the vortex wavelength s fflA aonpMMei d to 

Some new results for nonlinear vortices of ^TheSf and ^SriSent is shown to be 

experimental observations. The agreemen become important. For small wavelength 

vorticx.'the* nonli^ear^egime^of ^artlcul^ "the^ean^st^de ^which 

to“mL“r n“traf. Tlfe breaklown of this nonlinear state into a 
three-dimensional time dependent flow is also discussed. 


Journal of Fluid Mechanics. 

The interactions between longitudinal vOTtices^darcomj^yu^^ve^coi^dered^here 

axe strongly nonlinear, in the sense that e • ,, state N on li n ear interactions be- 

layer is completely SchSng waves are addressed first, and some analytical 

tween vortex flow and lollmien-bcnncnti g . i j oossibihty in the spatial- 

and computational propert.es are ^^-.J^cin^ displacement 

devdopmen. case toSiSf » —Mered for the com- 

thickness. Second, vortex/ Kayieig i cases Q f interest and some possible so- 

hi^n^jni^^Hes^Both^yp^^v^rtex/ToUmi^^ihchHng^^d.vortex/Raylei^h^ar^slKiH- 

l^feumbeS^nro^gly nonlinea'r nature is believed to make them very relevant 
to fully fledged transition to turbulence. 


Journal of Mechanics and Applied Mathematics. 

Recent theoretical work by Hall & ^ay«V« ^sSsSiWe 

high wavenumber Gortler vortices developing within a boundary layer now ar P 
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sec , on ^ry instability which takes the form of travelling waves confined to a thin region 
centered at the outer edge of the vortex. This work considered the case in which the secondary 
mode could be satisfactorily described by a linear stability theory and in the current paper our 
objective is to extend this investigation of Hall & Seddougui (1989) into the nonlinear Regime 

t n0t H*J y + t° e8 the SeCOndary . mode become nonlinear but it also 
Tnfitc Yl ,1 \ t0 modlfy the governing equations for the primary Gortler vortex 

In this case then, the vortex and the travelling wave drive each other and, indeed, the whole 

W^unde^T 18 de8cr ; b 5 d by a ? set of coupled, nonlinear differential equations. 

^®y eriake a Stuart- Watson type of weakly nonlinear analysis of these equations and 
m P arhcu _ lar - that on this basis there exist stable flow configurations in which the 
m ? m °, de I s of amplitude. Implications of our findings for practical situations are 

^ r ^ ed f and ,J t 18 S n°Z n that the fbeoretical conclusions drawn here are in good qualitative 
agreement with available experimental observations. 


TPA qp h R R ° UI ^ g M nd n S n S' S T ntharan ; U PP er semicontinuous global attractors for viscous flow. 
^thLatfcs ’ ’ ’ 28 Pag68 ‘ Submitted to Q^ly of Applied 


A particular form of regularization of the Navier-Stokes equations is studied. It has been 
shown that as the regularization parameter goes to zero, the global attractor for the regu- 
larized system converges to the global attractor for the conventional Navier-Stokes system 


HrfjSS 8 ’ J " u Uler in n A d ci a Z ier ~ Sto t 3 com P utations for two-dimensional geometries using 
unstructured meshes. ICASE Report No. 90-3 , January 2, 1990, 30 pages. Submitted to th? 
Canadian Aeronautics. 


A general purpose unstructured mesh solver for steady-state two-dimensional inviscid 
and viscous flows is described. The efficiency and accuracy of the method are enhanced 
y ^/-ultaneous use of adaptive meshing and an unstructured multigrid technique. A 
method for generating highly stretched triangulations in regions of viscous flow is outlined 

i^ n deflrSh°p C H edU R e f °u im P lei P entin S an algebraic turbulence model on unstructured meshes 
is described. Results are shown for external and internal inviscid flows and for turbulent 
viscous flow over a multi- element airfoil configuration. 


Nicol, David M. and Joel H. Saltz: An analysis of scatter decomposition. ICASE Report No 
2M, January 3, 1990, 20 pages. Submitted to IEEE Transactions on Computers. 1 


n -T^ 18 P a Per provides a formal analysis of a powerful technique known as scatter decom- 
position. Scatter decomposition divides an irregular computational domain into a large 

nrohahO 6qUal , S f ed pie ?f’ and distributes them modularly among processors. We use a 
probabilistic model of workload in one dimension to formally explain why, and when scatter 

d \^? S lh° n W °^ 8> • firSt re * SUlt i iS f hat if Correlation in workload is a convex function of 

d thCn SCatt n mg a m °f e fin ? ly decomposed domain yields a lower average processor 
workload variance. Our second result shows that if the workload process is stationary Gaus- 
sian and the correlation function decreases linearly in distance until becoming zero and then 
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remains zero, scattering a more finely decomposed domain yields a lower expected maximum 
processor workload. Finally we show that if the correlation function decreases linearly across 
the entire domain, then among all mappings that assign an equal number of domain pieces 
to each processor, scatter decomposition minimizes the average processor workload variance. 
The dependence of these results on the assumption of decreasing correlation is illustrated 
with situations where a coarser granularity actually achieves better load balance. 


Speziale, Charles G, Sutanu Sarkar, and Thomas B. Gatski: Modeling the pressure-strain 
correlation of turbulence - An invariant dynamical systems approach. lUAbb Keport INo_. 
9Q-5 . January 2, 1990, 52 pages. To be submitted to the Journal of Fluid Mechanics. 


The modeling of the pressure- strain correlation of turbulence is examined from a basic 
theoretical standpoint with a view toward developing improved second-order closure models. 
Invariance considerations along with elementary dynamical systems theory are used in the 
analysis of the standard hierarchy of closure models. In these commonly used models, the 
pressure-straincorrelation is assumed to be a linear function of the mean velocity gradien s 
with coefficients that depend algebraically on the anisotropy tensor. It is proven that for 
plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models 
is encapsulated by a relatively simple model which is only quadratically nonlinear in e 
anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the 
Launder, Reece, and Rodi model (as well as more recent models that have a considerably 
more complex nonlinear structure) in a variety of homogeneous turbulent flows. However, 
some deficiencies still remain for the description of rotating turbulent shear flows that are 
intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere 
introduction of more complex nonlmeanties. It is thus argued that the recent trend of 
adding substantially more complex nonlinear terms containing the anisotropy tensor may be 
of questionable value in the modeling of the pressure-strain correlation. Possible alternative 
approaches are discussed briefly. 


Iqbal M. Ashraf: Efficient algorithms for dilated mappings of binary trees. ICASE Report 
No. 90-6 . January 9, 1990, 41 pages. Submitted to IEEE Trans. Parallel & Distribute 
Computing. 


We address the problem of finding a 1-1 mapping of the vertices of a binary tree onto 
those of a target binary tree such that the son of a node on the first binary tree is mapped 
onto a descended of the image of the cost of this mapping, namely the dilation cost i.e. the 
maximum distance in the target binary tree between the images of vertices that are adjacent 
in the original tree. The other measure, expansion cost , is defined as the number ot extra 
nodes/edges to be added to the target binary tree in order to ensure a 1-1 mappmg. We 
describe an efficient algorithm to find a mapping of one binary tree onto another. We show 
that it is possible to minimize one cost of mapping at the expense of the other. 

This problem arises when designing pipelined Arithmetic Logic Units for special purpose 
computers. The pipeline is composed of ALU chips connected in the form of a binary tree. 
The operands to the pipeline can be supplied to the leaf nodes of the binary tree v ^ hl ^ th 
process and pass the results up to their parents. The final result is available at the root. As 
each new application may require a distinct nesting of operations, it is useful to be able to 
find a good mapping of a new binary tree over existing ALU tree. Another problem an 
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if every distinct required binary tree is known beforehand. Here it is useful to hardwire the 
pipeline in the form of a minimal supertree that contains all required binary trees. 


Koelbel, Charles, Piyush Mehrotra, and John Van Rosendale: Supporting shared data struc- 
tures on distributed memory architectures. ICASE Report No. 90-7 , January 9, 1990, 19 
pages. To appear in the Proceedings of the 2nd SIGPLAN Symposium on Principles and 
Practice of Parallel Programming, March 1990. 


Programming nonshared memory systems is more difficult than programming shared 
memory systems, since there is no support for shared data structures. Current programming 
languages for distributed memory architectures force the user to decompose all data struc- 
tures into separate pieces, with each piece “owned” by one of the processors in the machine, 
and with all communication explicitly specified by low-level message-passing primitives. This 
paper presents a new programming environment for distributed memory architectures, pro- 
viding a global name space and allowing direct access to remote parts of data values. We 
describe the analysis and program transformations required to implement this environment, 
and present the efficiency of the resulting code on the NCUBE/7 and IPSC/2 hypercubes. 


Jones, Mark T. and Merrell L. Patrick: Factoring symmetric indefinite matrices on high- 
performance architectures. ICASE Report No. 90-8 , January 9, 1990, 15 pages. Submitted 
to SIAM Journal for Matrix Analysis and Application. 


The Bunch-Kaufman algorithm is the method of choice for factoring symmetric indefi- 
nite matrices in many applications. However, the Bunch-Kaufman algorithm does not take 
advantage of high-performance architectures such as the Cray Y-MP. Three new algorithms, 
based on Bunch-Kaufman factorization, that take advantage of such architectures are de- 
scribed. Results from an implementation of the third algorithm are presented. 


Singer, Bart A., Thomas A. Zang, and Gordon Erlebacher: TS — Dean interactions in curved 
channel flow. ICASE Report No. 90-9 , January 18, 1990, 44 pages. Submitted to Journal 
of Fluid Mechanics. 


A weakly nonlinear theory is developed to study the interaction of TS waves and Dean 
vortices in curved channel flow. The predictions obtained from the theory agree well with 
results obtained from direct numerical simulations of curved channel flow, especially for 
low amplitude disturbances. At low Reynolds numbers the wave interaction is generally 
stabilizing to both disturbances, though as the Reynolds number increases, many linearly 
unstable TS waves are further destabilized by the presence of Dean vortices. 
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Durlofsky, Louis J. Stanley P»8«- Sub- 

hyperbolic conservation laws. ICASL Kepon 

milted to Journal of Computational Physics. 

A triangle based TVD 

tion of hyperbolic <mnservat.on laws m two space m a veraged data, which >s accom- 

the scheme lies in the nature of ^ tnt ° oolSion followed by a slope limiting procedure^ 
plished via a nearest neighbor linear in P resulting method is considerably more 

Two such Umiting procedures “ approximations which, like this scheme, ap- 

KT» P umerical result, for linear advection and 

Burgers’ equation are presented. 

Perthame, Benoit and ^^ASE Report 20 pag^ Submitted 

scalar conservation laws. 10 Aba tteporr i^j .» 

to SIAM Journal on Applied Mathematics. 

We construct a nonlinear 

DU a family of kinetic entropy ™ Equality on the macroscop, c 

ff-function, such that they .recov. =r K ^ov-^e > £ , he multidimensional case 

^ with e 1 0 to 

t" u^“?^cnof the corresponding conservation law. 

Optimization. 

One of the difficulties in designing controllers Galerkin 

attempting to calculate a state for the iP w m perform as designed when 

type approximations can be used to design co t No aasump tions, other than those 

on the approximating 

Osher, Stanley and Chi- Wang 

Hamilton- Jacobi equations. ICAbh rtepon 

mitted to SIAM Journal of Numerical Analysis. 

Hamilton- Jacobi (H-J) equations are to hyperbolic conserva- 
Ur fn o "dSSn fhe SSK simply the integrated version of the latter. 
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Sce ri s P ts: r t hi °i ° he,pM ." th = d -*” - «*- 

(ENO) schemes for H-J equations which yield uniform^ essentiall y non-oscillatory 
gions and resolve discontinuities in the in 8mooth re- 

procedure is adapted from that for hyperbolic conservation ItT schem ? construction 
schemes on a variety of one-dimensional and ■ , aw !\ We numerically test the 

related to control-optimization, and observe high orde^arr pro ^ ® ms> mcluding a problem 
resolution of discontinuities in the derives ^ *** 

cone.' ICASE Report* No .^9 0-1 4 ^ Fehr, . f^' hty ^"sonic flow past a sharp 
Comput ational Fluid Dynamic ? . ’ ’ °’ 6 ° pageS - Subm “ ed Theoretical l 

supersonic Sr^n““he« lSeriTTrvLture nl" 7 Y* Y* Y™ on a » har P «« in a 
This flow is then Analysed from the n n nT * Vr^ m the ph y sic . 8 of Problem. 
Indeed, the basic, non-axisymmetric disturbanrp »„ W ?• inear > temporal, inviscid stability, 
this class, and a so called “triply generalised” inflexio^ 10 ^^ 6 ^ lved , f ? r general flows of 
‘subsonic” neutral modes of instability. This C ? ndltl ? n 18 found for the existence of 

eralised inflexion condition found in planar flows a lth ** 18 ana !°8 ous to the well-known gen- 
depends on both axial and ariZtW “a" n2«1 8 ” PKSent case the “" di ‘ io " 

nuXr rf T 8 ~4Trfl SSlf” %l‘ aWily * at a *-*«m Mach 
instability may occur, peculiar toXws onSl tvoe -"Y 1 ? "T* 1 tbat a new mode of 
Additionally, asymptotic anX^abd ci^ to “Tf "T 

the cone are presented, and these give a Dartial ! P *- \ c< ? ne . / far downstream of 

mode of instability. g P rtla ( as y m Ptotic) description of this additional 


ulation ^ S an ?nln lS and Sim - 

Submitted to Theoretical and CompititSffi^i’ ^ ®’ 199 °’ 29 ^ 

to the general belief tha^suSTflows are^afcst^rom^ fc ^ r i bu l® nt M f ch numbe rs. Contrary 
velocity field remains small for all times) it is shown dlver g ence of the 7 

velocity field is negligibly small, it can grow rS ™ &t a- dlver « ence °f the initial 
is the inverse of the fluctuating Mach number & “°“' dl ™ en810 ual time scale which 

obtain a description of the flow in terms of its . w . hich enables one to 

has been developed to solve the initial value r>mh) 8 c 5^ ree an ^ vorticity-free components 
Mach number turbulent climes haw W cCswY', A s a ‘he various types of low 

Formulae are derived that accurately predict the level oT'co reapec ‘ mitral conditions, 
sients have disappeared. These reJSu ™ Verified^ u' ° f “ m P res J! bll l‘y after the initial tran- 
of isotropic turbulence. fied by ext ™"™ direct numerical simulations 
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Abarbanel Saul S Wai Sun Don, David Gottlieb, David H. Rudy, and James C. Townsend: 
Secondary' frequencies in the wake of a circular cylinder with vortex shedding. ICASE Repo rt 
No. 90-16 . February 9, 1990, 26 pages. Submitted to Journal of Fluid Mechanics. 

A detailed numerical study of two-dimensional flow past a circular cylinder at moder- 
ately low Reynolds numbers has been conducted using three different numerical algorithms 
fo r ^?he time-dependent compressible Navier-Stokes equations^ I was found -thrfrf 
the algorithm and associated boundary conditions were consistent and stable, tnen th 
or features of Se unsteady wake were well-predicted. However it was also found that even 
stable and consistent boundary conditions could introduce additional periodic phenomena 
reminiscent of the type seen in previous wind-tunnel experiments. However, these additional 
frequencies^ were ehnSnated by formulating the boundary conditions in terms of the charac- 
teristic variables 6 An analysis 7 based on a simplified model provides an explanation for this 

behavior. 


Decker Naomi H : On the parullel efficiency of the Frederickson-McBryan rnultigrid algo- 
rithm. ’ iCASE Report No. 90-17 , February 9, 1990, 10 pages. Submitted to SIAM Journa 

on Scientific and Statistical Computing. 


To take full advantage of the parallelism in a standard multigrid algorithm requires as 
methods. PSMG is simply a new and perhaps simpler was of achieving the same result 


to AIAA Journal. 

Theoretically based turbulence models have had success in predicting many features of 

that observed in experiments. Having thu, , «.nfirmed that 
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led to good agreement of the computed growth rates with the experimental data. Both the 
computations and the experiments indicate a dramatic reduction in the growth rate when 
the convective Mach number is increased. Experimental data on the normalized maximum 
turbulence intensities and shear stress also show a reduction with increasing Mach number. 
Ine computed values are in accord with this trend. 


Seddougui, S. 0., B °wles, and F. T. Smith: Surface - Cooling effects on compressible 
boundary-layer instability . ICASE Report No. 90-19 , February 28, 1990, 51 pages. Submit- 
ted to European Journal of Mechanics. 

,, Til f. influence of surface cooling on compressible boundary-layer instability is discussed 
theoretically for both viscous and inviscid modes, at high Reynolds numbers. The cooling 
enhances the surface heat transfer and shear stress, creating a high-heat-transfer sublayer 
Ihis has the effect of distorting and accentuating the viscous Tollmien-Schlichting modes to 
such and extent that their spatial growth rates become comparable with, and can even exceed , 
the growth rates of inviscid modes, including those found previously. This is for moderate 
cooling, and it applies at any Mach number. In addition, the moderate cooling destabilizes 
otherwise stable viscous or inviscid modes, in particular triggering outward-traveling waves 
at the of the boundary layer in the supersonic regime. Severe cooling is also discussed as 
it brings compressible dynamics directly into play within the viscous sublayer. All the new 
cooled modes found involve the heat-transfer sublayer quite actively, and they are often 
multi-structured in form and may be distinct from those observed in previous computational 
and experimental investigations. The corresponding nonlinear processes are also pointed out 
with regard to transition m the cooled compressible boundary layer. Finally, comparisons 
with Lysenko and Maslov s (1984) experiments on surface cooling are presented. 


Nicol, David M.: The cost of conservative synchronization in parallel discrete event simula- 
ICASE R eport No. 90-20 , May 9, 1990, 31 pages. Submitted to the Journal of the 

ACM. 


This paper analytically studies the performance of a synchronous conservative parallel 
discrete-event simulation protocol. The class of simulation models considered are oriented 
around a physical domain, and possess a limited ability to predict future behavior. Using a 
stochastic model we show that as the volume of simulation activity in the model increases 
relative to a fixed architecture, the complexity of the average per-event overhead due to 
synchronization, event list manipulation, lookahead calculations, and processor idle time 
approaches the complexity of the average per-event overhead of a serial simulation. The 
method is therefore within a constant factor of optimal. Our analysis demonstrates that 
on large problems-those for which parallel processing is ideally suited-there is often enough 
parallel workload so that processors are not usually idle. We also demonstrate the viability 
of the method empirically, showing how good performance is achieved on large problems 
using a thirty-two node Intel iPSC/2 distributed memory multiprocessor. 
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Nicol David M.: Performance bounds on 
port No. 90-21 , March 2, 1990, 29 pages, 
and Computer Simulations. 


parallel self-initiating discrete-event. ICASE 
Submitted to ACM Transactions on Modelling 


This paper considers the use of massively “HfafeSSg 

simulations of what we term self-mi la ing m . / D f other logical process, 

model schedules its own state re-evalua ion > P . .u re-evaluation. Our interest 
and sends its new state to other logical P^tmnkZn We co™d« the performance 
is in the effects of that communication on bounds on optimal perfor- 

of various synchronization protocols by d g PP t bounds on the performance 
mance, upper bounds on Time Warp’s performa^^ costs o{ 

of a new conservative protocol. Our ^alysis of Time WMP mci^.^ for ^ conservative 

state-saving and rollback. The ana P°' . . al>0 quan tifies the sensitivity of perfor- 

L r "°^ss^Coui>okS«d ability, and the probability distributions underlying 
the simulation. 


March 2.W9M5 pages. Submitted to IMA Journal of Numerical Analysis. 


The numerical simulation of steady planar 

pressible fluid through an abruptly contractmg channel us g P^^^^ fl 0 „ region 

methods is described. The key features of the method ar which are pointwise C> 

into a number of rectangukr subregfons and spectr ^ to the so i ut ion are obtained 

continuous across subregion interfaces. Spectral ^pp ^ ^ cQrner yortex decreases as 

for Reynolds numbers in the range [0, J* ^ 45 AS the Reynolds number is increased 
the Reynolds number increases fro detected downstream of the contraction at a 

“ th ' **"“* Mmber “ 

further. 
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Rosen, 


The convergence of solutions to the sXltons 

problem and associated Riccati equation *° r in , eaua tion as the length of the sampling 
to the corresponding continuous time problem ^^fis e^bhshed. Both the finite and 

interval (the sampling rate) tends °war ( . horizon case, strong continuity 

infinite time horizon problems are studied In the finite time ^honzon , & with a 

of the operators which define the control system an required. For the infinite 

stability and consistency “^Rion .on ei must be sta bilizable and detectable, 

STi.Cprf 

a flexible beam are presented and discussed. 
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thermod^™^ “T y for J ““. ideation of 

^^SSjS^£3£S£55SSS 

tlatHsg 

linew tSfne bTeH^ ^ *“? equat l ons ° f one dimensional linear thermoelasticity and a 
? pbne ^ ased 8 , ch j me f 18 J dlscuss ed and numerical results indicating how our aonroach 
might be used in a study of damping mechanisms in flexible structures tre presented h 

hiKhP^* 1 si^p^^ei^tr^^H^^ndn^^arbitraril^ 

g :4:KKSr. ir tr^No sssr and a signlficant speed up ° f 


Review of Fluid Mechanics, January 1991. 

„,±!> yll “ 1 .™ th ,° ds for tile development of Reynolds stress models in turbulence are 
d “ < j etMl - Zero > oue and two equation models are discussed along with second-order 
mZ.lfin “*■ CaS V’ ?; ade for , lhe su P erior P^ctive capabilities of siTnd-oX closure 

b^t^a^P^^^^'^^t^^^PS^Pcuu^a^d^inlmmogMwu^timlmlenc^ ^A^dLcussio^o^the 

“ Reynolds stras modeling is alao “ 
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Cai Wei David Gottlieb, and Ami Harten: Cell averaging Chebyshev methods for hyperbolic 
^bkmJ. ICASE Report No. 90-27 , March 28, 1990, 22 pages. To appear m Computer and 

Math with Application, special issue: Advance in P.D.E. 


This paper describes a cell averaging method for the Chebyshev approximations of first 
orderhyperbolic equations in conservation form. We present 

between pointwise data at the collocation points and cell averaged quantities, and vi 
versa This step trivial for the finite difference and Fourier methods is nontrivial for the 
global polynomials used in spectral methods. We then prove that the cell averaging .methods 
presented ^are stable for linear scalar hyperbolic equations and present num^ical simulations 
of shock-density wave interaction using the new cell averaging Chebyshev methods. 


Kami S.: Viscous shock profiles and primitive formulations. ICASE Report No. 90 : 28 

March 29, 1990, 24 pages. Submitted to SIAM Journal of Numerical Analysis. 


We consider weak solutions of hyperbolic systems in prinutive ( n °“: co 5 serva ^;nls OI Ibo f ck 
which a consistent conservation form exists. We show that for primitive formulat ons, shock 
mlat ion S C are ^ot 1 uniquely defined by the states to either side of the shock but also depend 
on the viscous path connecting the two. Scheme-dependent high order correction terms are 
derivedTrXce consistent viscous shock profiles. The resu ting primitive algorithm i 
conservative to the order of the approximation. One dimensional Euler calculations of flows 
containing strong shocks clearly show that conservation errors in primitive calculations are 
reduced to truncation levels and that both conservative and primitive flow calculations are 

of comparable quality. 
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ICASE INTERIM REPORTS 


Hi^andani, S^ma, /° e * Salt z , Harry Berryman, and Piyush Mehrotra: A scheme for sup- 
posing distributed data structures on multicomputers. Interim Report No. 9 January 4 1990 

Memory Computing Conference, ’ 


r ,^ e P ro P ose a data migration mechanism that allows an explicit and controlled mapping 
o data to memory. While read or write copies of each data element is assigned to a specific 
ocation in the memory of a particular processor. Our proposed integration of a data migra- 
tion scheme with a compiler is able to eliminate the migration of unneeded data that can 
occur m multiprocessor paging or caching. The overhead of adjudicating multiple concurrent 
writes to the same page or cache line is also eliminated. We present data that suggests that 
the scheme we suggest may be a practical method for efficiently supporting data migration. 
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ICASE COLLOQUIA 
October 1, 1989 through March 31, 1990 


Name/ Affiliation/Title D a ^ e 

Mr. Patrick Queutey, Laboratoire d’hydrodynamique Navale, France October 6 
“Three-Dimensional Incompressible Navier-Stokes 
Solutions over Bodies at Incidence” 

Professor J. Walker, Lehigh University October 13 

“Dynamics of Turbulent Wall Layers” 

Professor Max Gunzburger, Virginia Polytechnic Institute and State October 16 
University 

“Numerical Algorithms for the Boundary Control of Flows 
with Applications to Viscous Drag Minimization” 

Dr. Nahil A. Sobh, Old Dominion University October 18 

“Finite Element Computations on the Connection Machine” 

Dr. Tobias OrlofF, Institute for Computer Applications in Science October 19 
and Engineering 

“A Parallel 3D Triangle Renderer” 

Dr. Serge Petition, Yale University October 20 

“Subspace Methods as a Methodology to Solve Eigenproblems 
on Massively Parallel Architectures; Example of Arnoldi’s 
Method on the CM2” 

Professor Anthony Leonard, Rutgers University October 30 

“Dynamical Systems and Mixing in Fluid Mechanics” 

Professor Helen Reed, Arizona State University October 31 

“Receptivity of the Boundary Layer on a Semi- Infinite 
Flat Plate with an Elliptic Leading Edge” 

Professor Helen Reed, Arizona State University October 31 

“Stability of High-Speed Chemically Reacting 
and Three-Dimensional Boundary Layers” 

Dr. Andrew Neish, University College London November 6 

“Internal Hypersonic Nozzle Flow” 
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Name/ Affiliation/Title 


Date 


Dr. Gerald Hedstrom, Lawrence Livermore National Laboratory 
“The Behavior of Waves in Grids” 

Professor Swami Manohar, University of North Carolina, Charlotte 
“Large Parallel Architectures for Multigrid” 

Dr. Takashi Abe, The Institute of Space and Astronomical Science, 
Japan 

“Direct Simulation Monte Carlo Method in Rarefied 
Gas Flow Simulations” 

Dr. Naomi Decker, Institute for Computer Applications in Science 
and Engineering 

“On the Efficiency of the Frederickson-McBryan Multigrid Algorithm” 

Dr. John Goodrich, NASA Lewis Research Center 
“Hopf Bifurcation in the Driven Cavity” 

Professor Paul Reynolds, University of Virginia 
“Atomic Operations Without Locking” 

Professor Thomas F. Balsa, University of Arizona 
“The Evolution of Nonlinear Disturbances in a Supersonic 
Mixing Layer” 

Professor K. R. Rajagopal, University of Pittsburgh 
“Swirling Flows of Newtonian and Non-Newtonian Fluids 
and Some Related Problems” 

Dr. Abdelkader Frendi, Northeastern University 
“A Numerical Study of Ignition Limits for Methane- Air Mixtures” 

Professor Steven Orszag, Princeton University 
“Organization and Structure in Turbulent Flows” 

Mr. Kurt Bryan, University of Washington 
“An Inverse Problem in Impedance Tomography” 

Professor Bernard Grossman, VPI&SU 
“Anomalous Wave Speeds in Non-Equilibrium Flows 
with Multiple Translational Temperatures” 


November 8 


November 20 


November 28 


November 29 


December 5 


December 14 


December 18 


January 12 


January 19 
January 23 
February 1 


February 5 
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Professor George Karniadakis, Princeton University 
“Computation of Transitional and Turbulent Flows 


in Complex Geometries” 

Mr. Charles Koelbel, Purdue University 
“Compiling Programs for Nonshared Memory Machines” 


February 8 


Bracy Elton, University of California, Davis 
“Stability and Convergence Results for Lattice Gas„ 
and Lattice Boltzmann Methods for Solving PDEs 


February 9 


Mr. Steve Lee, University of Illinois at Urbana 
“Large Krylov Methods in the Solution of Nonlinear 
Initial Value Problems in ODEs and DAEs” 


February 12 


Professor Bernard Grossman, VPI&SU . . . „ 

“Integrated Aerodynamics/Structural Design Optimization 

Professor Otto Zeman, Stanford University ^ 

“Recent Developments in Modeling Compressible Turbulence 


February 13 


February 14 


Dr. Peter Jacobs, University of Queensland, Australia 
“Hypervelocity Flow in Axisymmetric Nozzles” 

Mr. Gregory Byrd, Stanford University 
“Streamline: Cache-Based Message Passing in Scalable 
Shared Memory Multiprocessors” 

Professor Biswa Datta, Northern Illinois University 
“Large-Scale and Parallel Computations in Control” 

Dr. Mark Morkovin, Professor Emeritus, Illinois Institute of 
Technology 

“On the Physics of Vorticity and Turbulence 
at Supersonic Speeds-Informal Ruminations” 

Professor W. 0. Criminale, University of Washington 
“The Evolution of Linearized Perturbations of Parallel Flows 


March 7 


March 8 


March 9 


March 12 


March 13 


Professor Thomas C. Corke, Illinois Institute of Technology 
“Three Dimensional Mode Resonance in Boundary Layers and Wa es 



Name/ Affiliation /Title 


Date 


Dr. Shahid H. Bokhari, University of Engineering & Technology, Lahore 
Pakistan ^ J 

“Pipelined Parallel Multigrid” 

Professor Joel Ferziger, Stanford University 

Direct Simulation of Turbulent Combusting Flows” 

Dr. Shahid H. Bokhari, University of Engineering & Technology, Lahore 
Pakistan * 

“A Network Flow Model for Analyzing Link Contention 
in Circuit-Switched Multicomputers” 


March 20 

March 21 
March 27 
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OTHER ACTIVITIES 


The Combustion Workshop co-sponsored by ICASE and NASA Langley Research Cen- 
ter was held October 2-4, 1989, at the Omni Hotel, Newport News, VA. Fifty-five peop e 
attended this workshop. The objectives of the workshop were to review the current status 
of knowledge on flame structure, flame stability, flame holding/extinction, ignition, chemi- 
cal kinetics, turbulence/kinetic interaction, transition to detonation, and reacting free shear 
layers. Invited speakers and their topics are listed below: 


Flame Structure: 

M. Smooke, Yale University 

B. Rogg, University of California, San Diego 

S. M. Correa, General Electric Corporate R and D Center 


Flame Stability: 

J. Buckmaster, University of Illinois 

T. Jackson, ICASE and Old Dominion University 

A. Kapila, Rensselaer Polytechnic Institute 

M. Matalon, Tech. Institute Engineering Science &: Applied Mathematics 

Mechanisms of Ignition, Extinction and Flame Stabilization 
in Subsonic and Supersonic Flows: 

C. K. Law, Princeton University 

A. Karagozian, University of California, Los Angeles 

I. Shepherd, Lawrence Berkeley Laboratory 

Mixing Power Concepts in Scramjet Combustor Design: 

J. R. Swithenbank, Sheffield University 
P. Givi, SUNY at Buffalo 

J. Riley, University of Washington 

Chemical Kinetics of Combustion Processes: 

R. W. Walker, Hull University 
F. L. Dryer, Princeton University 
R. Yetter, Princeton University 

C. J. Jachimowski, NASA Langley Research Center 

D. M. Golden, SRI International 
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The Interaction of Turbulence and Chemical Kinetics: 

W. Kollman, University of California, Davis 

D. R. Ballal, University of Dayton 

P. A. Libby, University of California, San Diego 

On the Transition from Deflagration to Detonation: 

J. Shepherd, Rennselaer Polytechnic Institute 
D. Pratt, University of Washington 
M. Sichel, University of Michigan 


A volume of the proceedings from this conference will be published by Springer-Verlag in 
the near future. 
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ICASE STAFF 


I. ADMINISTRATIVE 

Robert G. Voigt, Director 

Ph.D., Mathematics, University of Maryland, 1969 
Numerical Algorithms for Parallel Computers 

Linda T. Johnson, Office and Financial Administrator 

Etta M. Blair, Personnel/Bookkeeping Secretary 

Barbara A. Cardasis, Administrative Secretary 

Rosa H. Milby, Technical Publications/Summer Housing Secretary 

Shelly D. Millen, Technical Publications Secretary (Beginning December 1989) 

Cheryl A. Pruitt, Technical Publications Secretary (Through December 1989) 

Emily N. Todd, Executive Secretary /Visitor Coordinator 

II. SCIENCE COUNCIL for APPLIED MATHEMATICS and COMPUTER 
SCIENCE 

Tony Chen, Professor, Department of Mathematics, University of California at Los Angeles. 

John Hopcroft, Joseph C. Ford Professor of Computer Science, Cornell University. 

Anita Jones, Chairman, Department of Computer Science, University of Virginia. 

Robert MacCormack, Professor, Department of Aeronautics and Astronautics, Stanford Uni- 
versity. 

Joseph Oliger, Professor, Computer Science Department, Stanford University. 

Robert O’Malley, Jr., Chairman, Department of Mathematical Sciences, Rensselaer Poly- 
technic Institute. 

Stanley J. Osher, Professor, Mathematics Department, University of California. 

Werner C. Rheinboldt, Andrew W. Mellon Professor, Department of Mathematics and Statis- 
tics, University of Pittsburgh. 

John Rice, Chairman, Department of Computer Science, Purdue University. 

Burton Smith, Tera Computer Company, Seattle, WA. 
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Robert G. Voigt, Director, Institute for Computer Applications in Science and Engineering, 
NASA Langley Research Center. 

III. ASSOCIATE MEMBERS 

Saul S. Abarbanel, Professor, Department of Applied Mathematics, Tel-Aviv University, 
Israel. 

H. Thomas Banks, Professor, Center for Applied Mathematical Sciences, University of South- 
ern California. 

David Gottlieb, Professor, Division of Applied Mathematics, Brown University. 

Peter D. Lax, Professor, Courant Institute of Mathematical Sciences, New York University. 
Merrell L. Patrick, Professor, Department of Computer Science, Duke University. 


IV. CHIEF SCIENTIST 

M. Yousuff Hussaini - Ph.D., Mechanical Engineering, University of California, 1970. Com- 
putational Fluid Dynamics. (Beginning April 1978) 


V. LEAD COMPUTER SCIENTIST 

Joel H. Saltz - Ph.D., Computer Science, Duke University, 1985. Parallel Computing with 
Emphasis on Systems and Algorithmic Issues. (July 1989 to September 1992) 


V. SENIOR STAFF SCIENTIST 

Gordon Erlebacher - Ph.D., Plasma Physics, Columbia University, 1983. Computational 
Fluid Dynamics. (November 1989 to November 1994) 

Charles G. Speziale - Ph.D., Aerospace and Mechanical Sciences, Princeton University, 
1978. Fluid Dynamics with Emphasis on Turbulence Modeling and the Transition Process. 
(September 1987 to September 1991) 

John R. Van Rosendale - Ph.D., Computer Science, University of Illinois, 1980. Parallel 
Systems and Algorithms. (July 1989 to July 1991) 


VI. SCIENTIFIC STAFF 

Harry Scott Berryman - B.S., Computer Science, Memphis State, 1988. Performance Anal- 
ysis of Parallel Computing Systems. (July 1989 to July 1991) 

Richard G. Carter - Ph.D., Numerical Analysis, Rice University, 1986. Numerical Methods 
for Optimization Problems. (September 1987 to September 1990) 

Leon M. Clancy - B.S., Mechanical Engineering, University of Washington, 1971. System 
Manager. (December 1989 to Present) 
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Thomas W. Crockett - B.S., Mathematics, College of William and Mary, 1977. Parallel 
Systems Research. (February 1987 to February 1992) 

Naomi Decker - Ph.D., Mathematics, University of Wisconsin-Madison, 1987. Multigrid 
Method Applied to Hyperbolic Equations. (September 1987 to September 1990) 

Thomas M. Eidson - Ph.D., Mechanical Engineering, University of Michigan, 1982. Parallel 
Techniques for Computational Fluid Dynamics. (August 1989 to August 1991) 

Peter A. Jacobs - Ph.D., Computational Fluid Dynamics, University of Queensland, Aus- 
tralia, 1987. Compressible Fluid Dynamics. (January 1990 - January 1991) 

Fumio Kojima - Ph.D., Control Theory, Kyoto University, Japan, 1985. Probabilistic and 
Stochastic Methods for Optimal Control Problems. (September 1986 to September 1990) 

Dimitri Mavriplis - Ph.D., Mechanical and Aerospace Engineering, Princeton University, 
1988. Grid Techniques for Computational Fluid Dynamics. (February 1987 to February 
1992) 

Kirsten Morris - Ph.D., Electrical Engineering, University of Waterloo, Canada, 1989. Con- 
trol Theory. (September 1989 to September 1991) 

Tobias B. Orloff - Ph.D., Mathematics, Princeton University, 1983. High Performance 
Graphics. (September 1988 to October 1989) 

Yuh-Roung Ou - Ph.D., Aerospace Engineering, University of Southern California, 1988. 
Control Systems for Fluid Dynamics. (November 1988 to November 1990) 

Douglas H. Peterson, Jr. - M.S., Computer Science, College of William and Mary, 1987. 
Computer Systems and Networking. (October 1985 to December 1989) 

Peter W. Protzel - Ph.D., Electrical Engineering, Technical University of Braunschweig, 
Germany, 1987. Reliability of Computing Systems. (March 1987 to September 1990) 

Sutanu Sarkar - Ph.D., Mechanical and Aerospace Engineering, Cornell University, 1988. 
Fluid Dynamics, Turbulence Modeling, Compressible Turbulence. (September 1988 to Septem- 
ber 1990) 

Jeffrey S. Scroggs - Ph.D., Computer Science, University of Illinois at Urbana, 1988. Domain 
Decomposition for Differential Equations. (July 1988 to July 1990) 

Sharon 0. Seddougui - Ph.D., Applied Mathematics, University of Exeter, England, 1988. 
Compressible Fluid Dynamics. (June 1988 to June 1990) 


VII. VISITING SCIENTISTS 

Takashi Abe - Ph.D., Aeronautics, University of Tokyo, 1978. Associate Professor, Institute 
of Space and Astronautical Science, Kanagawa, Japan. Monte Carlo Methods. (November 
to December 1989) 

Shahid Bokhari - Ph.D., Electrical and Computer Engineering, University of Massachusetts, 
Amherst, 1978. Professor, Electrical Engineering, University of Engineering & Technology, 
Lahore, Pakistan. Parallel Computing Systems. (October 1989 - January 1991) 
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Dave Broutman - Ph.D., Physical Oceanography, Scripps Institute Oceanography/University 
of California at San Diego, 1982. Lecturer/Research Fellow, Department of Mathemat- 
ics, University of New South Wales, Sydney, Australia. Multi-domain Spectral Methods. 
(December 1989 to January 1990) 

Peter W. Duck - Ph.D., Fluid Mechanics, Univeristy of Southampton, United Kingdom, 
1975. Lecturer in Mathematics, Department of Mathematics, University of Manchester, 
United Kingdom. Numerical Solution of Unsteady Boundary Layer Equations. (January 
1990) 

Michael Gaster - Ph.D., Aerodynamics, London University, 1963. Professor of Aeronautical 
Engineering, Department of Engineering, Cambridge University. Instability and Transition. 
(February to March 1990) 

Bernard Grossman - Ph.D., Astronautics, Polytechnic Institute of Brooklyn, 1969. Professor 
and Associate Department Head, Department of Aerospace and Ocean Engineering, Virginia 
Polytechnic Institute and State University. Multidisciplinary Design. (January to March 
1990) 

Philip Hall - Ph.D., Mathematics, Imperial College, England, 1973. Professor, Department 
of Applied Mathematics, University of Exeter, England. Computational Fluid Dynamics. 
(November to December 1989 and January to February 1990) 

Piyush Mehrotra - Ph.D., Computer Science, University of Virginia, 1982. Assistant Pro- 
fessor, Department of Computer Science, Purdue University. Programming Languages for 
Multiprocessor Systems. (January to September 1990) 

Serge G. Petiton - Ph.D., Computer Science and Applied Mathematics, Pierre and Marie 
Curie University (Paris VI), 1988. Post Doctoral Fellow, Research Center for Scientific 
Computation, Yale University. Parallel Methods for Eigenvalue Problems for Non-Symmetric 
Linear Systems. (October 1989) 

Eitan Tadmor - Ph.D., Numerical Analysis, Tel-Aviv University, 1979. Senior Lecturer, 
Department of Applied Mathematics, Tel-Aviv University. Numerical Methods for Partial 
Differential Equations. (July 1989 - July 1990) 

VIII. CONSULTANTS 

Loyce M. Adams - Ph.D., Applied Mathematics, University of Virginia, 1983. Associate Pro- 
fessor, Department of Applied Mathematics, University of Washington. Numerical Methods 
for Parallel Computing Systems. 

Alvin Bayliss - Ph.D., Mathematics, New York University, 1975. Associate Professor, Tech- 
nological Institute, Northwestern University. Numerical Methods for Partial Differential 
Equations. 

Marsha J. Berger - Ph.D., Numerical Analysis, Stanford University, 1982. Research Asso- 
ciate, Courant Institute of Mathematical Sciences. Numerical Methods for Partial Differen- 
tial Equations. 

Percy Bobbitt - B.S., Aeronautics, Catholic University of America, 1949. NASA Langley 
Research Center - Retired. Fluid Mechanics. 
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Achi Brandt - Ph.D., Mathematics, Weizmann Institute of Science, 1965. Professor, Applied 
Mathematics Department, Weizmann Institute of Science, Israel. Multigrid Methods. 

Dennis W. Brewer - Ph.D., Mathematics, University of Wisconsin, Madison, 1975. Asso- 
ciate Professor, Department of Mathematical Sciences, University of Arkansas. Methods for 
Parameter Identification and Estimation. 

John A. Burns - Ph.D., Mathematics, University of Oklahoma, 1973. Professor, Virginia 
Polytechnic Institute and State University. Numerical Methods in Feedback Control and 
Parameter Estimation. 

Peter R. Eiseman - Ph.D., Mathematics, University of Illinois, 1970. Senior Research Sci- 
entist and Adjunct Professor, Department of Applied Physics and of Nuclear Engineering, 
Columbia University. Computational Fluid Dynamics. 

Robert E. Fennell - Ph.D., Mathematics, University of Iowa, 1969. Professor, Department 
of Mathematical Sciences, Clemson University. Control Theory for Multivariable Systems. 

James F. Geer - Ph.D., Applied Mathematics, New York University, 1967. Professor, Sys- 
tems Science and Mathematical Sciences, Watson School of Engineering, Applied Science 
and Technology, SUNY-Binghamton. Perturbation Methods and Asymptotic Expansions of 
Solutions to Partied Differential Equations. 

J. Steven Gibson - Ph.D., Engineering Mechanics, University of Texas at Austin, 1975. Asso- 
ciate Professor, Department of Mechanical, Aerospace and Nuclear Engineering, University 
of California at Los Angeles. Control of Distributed Systems. 

Chester E. Grosch - Ph.D., Physics - Fluid Dynamics, Stevens Institute of Technology, 
1967. Professor, Department of Computer Science and Slover Professor, Department of 
Oceanography, Old Dominion University. Hydrodynamic Stability, Computational Fluid 
Dynamics, Unsteady Boundary Layers and Algorithms for Array Processors. 

Philip Hall - Ph.D., Mathematics, Imperial College, England, 1973. Professor, Department 
of Applied Mathematics, University of Exeter, England. Computational Fluid Dynamics. 

Amiram Harten - Ph.D., Mathematics, New York University, 1974. Associate Professor, 
Department of Mathematics, Tel-Aviv University, Israel. Numerical Solution for Partial 
Differential Equations. 

Thorwald Herbert - Ph.D., Aerospace Engineering, University of Stuttgart, Germany, 1978. 
Professor, Department of Mechanical Engineering, Ohio State University. Fluid Dynamics. 

Kazufumi Ito - Ph.D., Systems Science and Mathematics, Washington University, 1981. 
Assistant Professor, Department of Mathematics, University of Southern California. Control 
Theory. 

Thomas L. Jackson - Ph.D., Mathematics, Rensselaer Polytechnic Institute, 1985. Assistant 
Professor, Department of Mathematics, Old Dominion University. Numerical and Analytical 
Methods for Chemically Reacting Flows. 

Antony Jameson - Ph.D., Magnetohydro-Dynamics, Cambridge University, England, 1963. 
James S. McDonnell Distinguished Professor, Department of Mechanical and Aerospace 
Engineering, Princeton University. Computational Fluid Dynamics. 
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Charles R. Johnson - Ph.D., Mathematics and Economics, California Institute of Technology, 
1972. Professor, Department of Mathematics, College of William and Mary. Numerical 
Linear Algebra. 

Mark T. Jones - B.S., Computer Science, Clemson University, 1986. Teaching Assistant, 
Duke University. Parallel Algorithms for Numerical Linear Algebra. 

Harry F. Jordan - Ph.D., Physics, University of Illinois, 1977. Professor, Department of 
Electrical and Computer Engineering, University of Colorado at Boulder. Parallel Compu- 
tation. 

Ashwani K. Kapila - Ph.D., Theoretical and Applied Mechanics, Cornell University, 1975. 
Associate Professor, Department of Mathematical Sciences, Rensselaer Polytechnic Institute. 
Ordinary and Partial Differential Equations, Asymptotic Methods. 

Stephen L. Keeling - Ph.D., Mathematics, University of Tennessee, 1986. Assistant Professor, 
Vanderbilt University. Analysis of Mathematical Models for Problems Arising in Acoustics. 

Edward J. Kerschen - Ph.D., Mechanical Engineering, Stanford University, 1978. Associate 
Professor, Department of Aerospace and Mechanical Engineering, University of Arizona. 
Flow Dynamics. 

David E. Keyes - Ph.D., Applied Mathematics, Harvard University, 1984. Assistant Pro- 
fessor, Mechanical Engineering, Yale University. Parallelizaton of Numerical Procedures 
Appropriate for the Study of Combustion. 

Heinz-Otto Kreiss - Ph.D., Mathematics, Royal Institute of Technology, Sweden, 1960. Pro- 
fessor, Department of Applied Mathematics, California Institute of Technology. Numerical 
Analysis. 

William D. Lakin - Ph.D., Applied Mathematics, University of Chicago, 1968. Professor, 
Department of Mathematics and Statistics, University of Vermont. Computational Fluid 
Dynamics. 

David G. Lasseigne - Ph.D., Applied Mathematics, Northwestern University, 1985. Assistant 
Professor, Department of Mathematics and Statistics, Old Dominion University. Computa- 
tional Fluid Dynamics. 

Anthony Leonard - Ph.D., Nuclear Engineering, Stanford University, 1963. Professor of 
Aeronautics, California Institute of Technology. Fluid Physics. 

Randall J. LeVeque- Ph.D., Computer Science, Stanford University, 1982. Assistant Profes- 
sor Department of Mathematics, University of Washington. Numerical Solution of Partial 
Differential Equations. 

John L. Lumley - Ph.D., Aeronautics, John Hopkins University, 1957. Professor, Department 
of Mechanical and Aerospace Engineering, Cornell University. Mathematical Aspects of 
Turbulence. 

Robert W. MacCormack - M.S., Mathematics, Stanford University. Professor, Department 
of Aeronautics and Astronautics, Stanford University. Computational Fluid Dynamics and 
Numerical Analysis. 

Ravi Mirchandaney - Ph.D., Computer Engineering, University of Massachusetts, 1987. Re- 
search Scientist, Department of Computer Science, Yale University. Parallel Run-time Sup- 
port Systems. 
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Mark V. Morkovin -Ph.D., Applied Mathematics, University of Wisconsin, 1942. Profes- 
sor Emeritus, Department of Mechanical and Aerospace Engineering, Illinois Institute of 
Technology. Transition Process in Aerodynamics. 

Katherine A. Murphy - Ph.D., Applied Mathematics, Brown University, 1983. Assistant 
Professor, Department of Mathematics, University of North Carolina, Chapel Hill. Control 
Theory and Estimation of Parameters. 

David M. Nicol - Ph.D., Computer Science, University of Virginia, 1985. Professor, Depart- 
ment of Computer Science, College of William and Mary. Mapping Algorithms onto Parallel 
Computing Systems. 

Roy A. Nicolaides - Ph.D., Mathematics, University of London, 1972. Professor, Depart- 
ment of Mathematics, Carnegie- Mellon University. Numerical Solution of Partial Differential 
Equations. 

James M. Ortega - Ph.D., Mathematics, Stanford University, 1962. Professor and Chairman, 
Department of Applied Mathematics, University of Virginia. Numerical Methods for Partial 
Differential Equations. 

Stanley J. Osher - Ph.D., Functional Analysis, New York University, 1966. Professor, Depart- 
ment of Mathematics, University of California at Los Angeles. Methods for the Numerical 
Analysis of Partial Differential Equations. 

Ugo Piomelli - Ph.D., Mechanical Engineering, Stanford University, 1987. Professor, Depart- 
ment of Mechanical Engineering, University of Maryland. Subgrid Scale Reynold’s Stress 
Modelling and Large Eddy Simulation of Turbulent Flows. 

Terrence W. Pratt - Ph.D., Mathematics/Computer Science, University of Texas at Austin, 
1965. Professor, Department of Computer Science, University of Virginia. Programming 
Languages. 

Daniel A. Reed - Ph.D., Computer Science, Purdue University, 1983. Assistant Professor, 
Department of Computer Science, University of Illinois. Parallel Processing. 

Helen L. Reed - Ph.D., Engineering Mechanics, Virginia Polytechnic Institute and State 
University, 1981. Associate Professor, Department of Mechanical Engineering, Arizona State 
University. Computional Fluid Dynamics. 

Eli Reshotko - Ph.D., Aeronautics and Physics, California Institute of Technology, 1960. In- 
terim Dean, Case Western Reserve University. High Speed Aerodynamics with an Emphasis 
on Transition, Turbulence and Combustion. 

Paul F. Reynolds - Ph.D., Computer Science, University of Texas at Austin, 1979. Assistant 
Professor, Department of Computer Science, The University of Virginia. Parallel Computing 
Systems. 

I. Gary Rosen - Ph.D., Applied Mathematics, Brown University, 1980. Assistant Professor, 
Department of Applied Mathematics, University of Southern California. Numerical Methods 
for Problems in Control Systems. 

Paul E. Saylor - Ph.D., Mathematics, Rice University, 1986. Associate Professor, Computer 
Science Department, University of Illinois, Urbana. Iterative Solution of Linear Algebraic 
Equations and Algorithms for Parallel Computers. 
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Chi-Wang Shu - Ph.D., Mathematics, University of California, Los Angeles, 1986. Assistant 
Professor, Division of Applied Mathematics, Brown University. Partial Differential Equa- 
tions. 

Lawrence Sirovich - Ph.D., Fluid Mechanics, John Hopkins Univesity, 1960. Professor, Di- 
vision of Applied Mathematics, Brown University. Fluid Mechanics. 

Edward A. Spiegel - Ph.D., Astronomy, University of Michigan, 1958. Rutherford Professor 
of Astronomy, Department of Astronomy and Astrophysics, Columbia University. Mathe- 
matical Aspects of Turbulence. 

Katepalli R. Sreenivason - Ph.D., Aeronautical Engineering, Indian Institute of Science, 
1975. Professor and Chairman, Department of Mechanical Engineering, Yale University. 
Transition and Turbulence. 

Shlomo Ta’asan - Ph.D., Applied Mathematics, Weizmann Institute, 1984. Scientist, De- 
partment of Applied Mathematics, The Weizmann Institute of Science, Israel. Multigrid 
Methods for Partial Differential Equations. 

Saleh Tanveer - Ph.D., Applied Mathematics, California Institute of Technology, 1984. Pro- 
fessor, Department of Mathematics, Ohio State University. Problems for Crystal Growth. 

Lu Ting - Ph.D., Aeronautics, New York University, 1951. Professor, Courant Institute of 
Mathematical Sciences, New York University.. Fluid Mechanics. 

Sherryl J. Tomboulian - Ph.D., Computer Science, Duke University, 1986. Maspar Computer 
Corporation, Sunnyvale, California. Parallel Computing Systems Design. 

Eli Turkel - Ph.D., Applied Mathematics, New York University, 1970. Associate Profes- 
sor, Department of Applied Mathematics, Tel-Aviv University, Israel. Computational Fluid 
Dynamics. 

Bram van Leer - Ph.D., Theoretical Astrophysic, Leiden State University, The Netherlands, 
1970. Professor, Department of Aerospace Engineering, University of Michigan. Computa- 
tional Fluid Dynamics. 

Yun Wang - M.S., Applied Mathematics, University of Southern California, 1986. Research 
Associate, Center for Control Sciences, Division of Applied Mathematics, Brown University. 
Control Theory. 

J. Christian Wild - Ph.D., Computer Science, Rutgers University, 1977. Assistant Profes- 
sor, Department of Computer Science, Old Dominion University. Concurrent Computing 
Systems. 

IX. STUDENT ASSISTANTS 

Michael Arras - Graduate Student at College of William and Mary. (October 1988 to Present) 

X. GRADUATE FELLOWS 

Charles H. Koelbel - Graduate Student at Purdue University. (March to May 1990) 

Haim Nessyahu - Graduate Student at Tel-Aviv University, Israel. (October to December 
1989) 

Tamir Tassa - Graduate Student at Purdue University. (October to December 1989) 
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